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THE MANUAL 


PLOT 50: STATISTICS, VOLUME 1 is a collection of general purpose programs categorized 
as Distributions, Tabled Values, Tests and Simple Regressions, and One Sample Analysis. 
The programs have been designed to perform rapid calculations and to provide accurate re- 
sults in easily usable form. This manual is a guide to the execution of the programs. 


Although not required to execute programs described here, the 4050 series manuals listed 
below provide much valuable information: 


—À 


TITLE TEKTRONIX PART NUMBER 
4051 Graphic System Operator's Manual 070-1940-00 
n 4051 Graphic System Reference Manual 070-2056-00 
Plot 50: Introduction To Programming 070-2058-00 
In BASIC 
q Plot 50: Introduction To Graphics 070-2059-00 


Programming In BASIC 


) 


Each program is treated as a separate entity, and information applicable to more than one 
program is given for each program to which it applies. This repetition eliminates the need 

for “page-flipping”, and any program can be executed if you follow the instructions con- 

tained here without referring to other sources. — 


pão | 


An introduction, a list of assigned variable names, operating instructions, and examples are 
provided for each program. The introduction contains a description of the program, input 
requirements, expected output, hardware requirements, and program limits. Methods used 
to derive results are explained briefly, and supporting references are listed. 


Lists of variables used in each program are provided but are not needed for program 
execution. In the lists, the word “scratch” is used to define interim working values which 
are used in program operations but are not displayed as output. The variable lists may be 
disregarded unless program modification is undertaken. 


fas d en 


a | 


Program loading and execution are described in the OPERATING INSTRUCTIONS which 
are brief, explicit, and in logical sequence. Throughout the instructions, messages which are 
displayed by the system are shown in upper case letters. Upper case is also used to designate 
specific keys on the keyboard (e.g., RETURN, BREAK, AUTO LOAD, etc.). 
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USER-DEFINABLE KEYS 


There are ten user-definable keys at the upper left corner of the 4051 keyboard. Each key 
may have a specific value when pressed alone, and another value when pressed with the 
SHIFT key down. When a program which uses the keys has been read into memory, the 
values assigned to the keys are associated with designated program functions. Then, if a 
user-definable key is pressed, with our without the SHIFT key, the corresponding program 
function is activated. 


LI ELS O E 


Overlay templates for the user-definable keys are included with the STATISTICS, VOLUME 
1 package, and information concerning their use is part of all applicable program execution 
instructions. 


ES ES ES * 


OPERATING CONSIDERATIONS 


During program execution, a question mark flashing on the screen indicates that the program 
is waiting for input data values. 


EC. 


When the display screen is full, a flashing “F” appears in the upper left corner. No additional į y 
data will be displayed or accepted until the screen is cleared by pressing the PAGE key. esf 


Messages which appear because of error conditions may be avoided by following the operating i 
instructions carefully. An explanation of the messages is contained in the Graphic System 
Reference manual. T 


Whenever possible, the programs in STATISTICS, VOLUME 1 avoid rigid syntactical rules, 
and human error is forgiven or ignored when practical. This flexibility allows minor variation 
in program execution procedure. Therefore, no claim is made that the methods detailed 
here are the "only" way. The instructions will lead to effective results however, and they 

are offered as a sound base for development of individualized techniques. 
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F-DISTRIBUTION 


DESCRIPTION 

This program will calculate the significance level (right-tail probability) of a designated F- 
value. Five digit accuracy wil! be determined for significance levels between .0001 and .9999 
when degrees of freedom are 1000 or less. Accuracy of two digits or more will result when 
both degrees of freedom are greater than 1000. 

A positive F-value and degrees of freedom for the numerator and denominator are required 
input. Right-tail probability is the output. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


The F-value must be positive. 


Degrees of freedom may be any real number greater than 0.01. 


METHODS 


The density function of the F-distribution with n and d degrees of freedom is: 


n^ /2qd/24ln -2)/2 (gt) -1/21n*d), 
f(n,d,t) = ———— 
B(n/2,d/2) 


. t» 0 where B(n/2,d/2) is the complete beta. 
The calculated right tail probability is S3 f(n,d,t) dt. 
The F-value is transformed to a beta value. 
When both degrees of freedom are 1000 or less, a continued fraction is used to evaluate 


significance levels (References 1, 2). If both are greater than 1000, a normal approximation 
is used (Reference 3). 
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F-DISTRIBUTION I 
REFERENCES L 
1. Aroian, T. A., "Continued Fractions for the Incomplete Beta Function”, Annals of 
Mathematical Statistics, 12 (1941), pp 218—223. L 
2. Kopitzke, R. W., “Unpublished Notes”. 
3. Peizer, D. B., Pratt, J. W., "A Normal Approximation for Binomial, F, Beta, and Other LU 
Common Related Tail Probabilities, | “Journal of the American Statistical Association, 
63 (1968), pp.1457— 1483. U 


VARIABLES 


1-2 


VARIABLE 


Denominator degrees of freedom/2 
Scratch 

Scratch ; 
Numerator degrees of freedom/2 V | 
Scratch 

Scratch 

Value of F 

Numerator degrees of freedom 


(— 


Denominator degrees of freedom 
Scratch 

Value of beta 

Scratch 

Scratch 


Scratch E 
Scratch 


Scratch L 


E 
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PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the key- 
board. 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message 
ENTER THE PROGRAM NUMBER YOU WANT. 


3. Enter 1; press RETURN. 
4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 2; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the I/O light goes out, execute the program. 


Rees 


PROGRAM EXECUTION 


1. NUM.DF, DEN.DF, F is displayed on the screen. Enter the degrees of freedom and the 
F-value, using commas to separate values. Press RETURN. 


2. The significance level is displayed. 


3. Other values may be entered by repeating step 1. To exit from the program, press the 
RETURN key twice. 


nes ood 


EXAMPLE 


=3 


F DISTRIBUTION 


NUM. DF, DEN. DF, F 5,15,2.4 
NUM. DF =5 DEN. DF =15 P(F)=2.4) = 0.09660? 


NUM. DF, DEN. DF, F 158,600.1.2 
NUM. DF =150 DEN. DF -600 P(F>=1.2> = 0.07191 


NUM. DF, DEM. DF, F 2,2,19 
NUM. DF =2 DEN. DF =2 PcF>=19) = 0.05 


^N NUM. DF, DEN. DF, F 50,20,.74 
) NUM. DF =50 DEN. DF =28 P(F>=0.74) = 8.89729 


1 NUM. DF, DEN. DF, F 
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DESCRIPTION 


This program calculates significance levels for a one-tail and two-tail test for the Students 
t-value. Accuracy of five or more digits will result for significance levels between .0001 and 
.9999 when degrees of freedom are 1000 or less. For degrees of freedom greater than 1000, 


accuracy will be at least two digits. 


The degrees of freedom and a calculated Students t-value are required input. Input values 


are displayed with the resultant significance levels. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


Degrees of freedom should be a real number greater than 0.5. 


METHODS 


The density function for the t-distribution with n degrees of freedom is: 


fina) = T ((n+1)/2) * (1+x2 /n) 7 (nt 1)/2 


VIInP(n/2) 


where: 


—oo << x< 00 + 


For degrees of freedom n and calculated t-value t, the quantities calculated are given by: 


Letr = It]: 


The two-tail significance is: 


1—s" fin,x) dx * 
—r 


The one-tail significance is: 


I f(n,x) dx = Pm f(n,x) dx * 


The calculated t-value is transformed to a beta variable. 
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If degrees of freedom are 1000 or less, a continued fraction is used to evaluate significance 
levels (References 1 and 2). If degrees of freedom are greater than 1000, a normal approxi- 
mation is used (Reference 3). 


onal 


REFERENCES 


1. Aroian, T. A., "Continued Fractions for the Incomplete Beta Function", Annals of 
Mathematical Statistics, 12 (1941), pp 218—223. 

2. Kopitzke, R. W., "Unpublished Notes". 

3. Peizer, D. B., Pratt, J. W., "A Normal Approximation for Binomial, F, Beta, and Other 


Common Related Tail Probabilities, |" Journal of the American Statistical Association, 
63 (1968), pp 1457—1483. 
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VARIABLES 


) 


VARIABLE USAGE 


Degrees of freedom/2 


t A1 Scratch 
! A2 Scratch 
B Degrees of freedom/2 
n B1 Scratch 
B2 Scratch 
n F Value of t 
: F2 Scratch 
Counter 


3 
SE 


Value of beta 


T Scratch 
1 U Scratch 
l V1 Scratch 
X Scratch 
Y4 Scratch 
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PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY U 
1. Insert the program tape. U 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 2; press RETURN. i 
4. When the 1/0 light goes out, execute the program. 


EI 


MANUAL LOADING FROM THE KEYBOARD 


1. Insert the program tape. i 
2. Enter FIND 3; press RETURN. 

| 
3. Enter OLD; press RETURN. VW 1 


4. When the I/O light goes out, execute the program. 


PROGRAM EXECUTION i 
1. When DEGREES OF FREEDOM, T is displayed, enter the values (separated by a 

comma) and press RETURN. E 

2. One-tail and two-tail significance levels are displayed. i 


3. Repeat the operation by entering other values, or exit from the program by pressing ; 
the BREAK key twice. I 


EXAMPLE i 
T DISTRIBUTION 


DEGREES OF FREEDOM ,T 65,1.4 
N=65 T=1.4 ONE TAIL=0.083134 TWO TAIL=0. 166268 


DEGREES OF FREEDOM ,T 3,-1.5 
N25 T2-1.5 ONE TAIL=0.096952 TWO TAIL=6. 193984 


pi OF FREEDOM «T ee 1.96 
1200 T=1.96 ONE TAIL=0.025113 TWO TAIL=B. 050226 


hegg OF FREEDOM sT 690,1.7 L 


Ts1.7 ONE TAIL=0.044731 TWO TAIL=0. 089582 
DEGREES OF FREEDOM ,T NE 


) 


D 29 4) 3 
p 


— 


ved 2 032 


DISTRIBUTIONS 


CHI-SQUARE DISTRIBUTION 


DESCRIPTION 


The chi-square distribution program calculates the right-tail probability (significance level) 
of the chi-square distribution. For degrees of freedom greater than 0.5 and significance 
levels between .0001 and .9999, probabilities are calculated with at least five figures of 
accuracy. 


Degrees of freedom and a chi-square value are required input. These values and the signifi- 


cance level are the output. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


Degrees of freedom may be any real number greater than 0.5. 


METHODS 


The density function of a chi-square random variable with n degrees of freedom is: 


fint) = et/2 ¢n/2 —1 . 


1 
P(n/2) 20/2 
The significance level calculated for a chi-square value x with n degrees of freedom is: 


P(X>x) = Pe f(n,t) dt 


If the degrees of freedom are 200 or less, two different continued fractions are used (Refer- 
ences 1 and 3). If the degrees of freedom are greater than 200, a normal approximation is 
used (Reference 2). 
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REFERENCES 


1. Khovanskii, A. N., The Applications of Continued Fractions and Their Generalizations 
to Problems in Approximation Theory, translated by P. Wynn (1963), Groningen: 


P. Noordhoff. 


2. Peizer, D. B., Pratt, J. W., "A Normal Approximation for Binomial, F, Beta, and Other 


Common Related Tail Probabilities, |, Journal of the American Statistical Association, 
63 (1968), pp 1416—1456. 


3. Wall, H. S., Analytic Theory of Continued Fractions, New York: D. Van Nostrand, 1948. 


VARIABLES 


VARIABLE 


PROGRAM LOADING 


USAGE 


Degrees of freedom/2 
Scratch 
Scratch 
Chi-square/2 
Scratch 
Scratch 


Scratch 

Chi-square 

Counter 

Degrees of freedom 
Scratch 

Scratch 

Scratch 


This program may be loaded automatically from the directory, or manually from the 


keyboard. 
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CHI-SQUARE DISTRIBUTION 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 3; press RETURN. 
4. When the I/O light goes out, execute the program. 


-3 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 4; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the I/O light goes out, execute the program. 


-J2 = 


-3 


/ ^ PROGRAM EXECUTION 


1. When the displayed message is DEGREES OF FREEDOM, CHI-SQUARE, enter the 
values and press RETURN. The values should be separated by a comma. 


2. The input values and the calculated significance level are displayed. 


Es 


3. Other values may be entered for further calculations, or the BREAK key may be pressed 
twice to exit from the program. 


ssa 252 


EXAMPLE 
1 
CHI-SQUARE DISTRIBUTION 
n DEGREES OF FREEDOM,CHI-SQUARE 50,62.3 
| Ha 50 PCCHI SQUARE 262.3 > «0.11372 
DEGREES OF FREEDOM, CHI-SQUARE 12,49,5 
1 12 — PKCHI SQUARE 549.5 > =1.7099E-6 
| GREES OF FREEDOM,CHI-SQUARE 12,2.3 
DEGREES OPE CHI SQUARE 92.2 > «0.99879 
! DEGREES OF FREEDOM, CHI-SQUARE 1500, 1472 
AN Ne 1580 PCCHI SQUARE >1472 ) =0.69225 
1 | DEGREES OF FREEDOM, CHI-SQUARE 
n Plot 50: Stat Vol 1 l 1 -9 
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GAMMA DISTRIBUTION L 
DESCRIPTION i 
This program calculates the right tail probability for a two parameter gamma distribution. 

When significance levels are between .0001 and .9999, probability statements are calculated 

with a minimum 3 digit accuracy for all values of a/b greater than .25. U 

Parameters a and b and the gamma value are required input. || 

HARDWARE REQUIREMENTS L 

A Tektronix 4050-Series Graphic System is required. 

PROGRAM LIMITS L 

The ratio of a to b must be greater than .25 (a/b > 1/4). i 

METHODS uw 

The density function of the gamma distribution with parameters a and b is: i as 
1 gx L 

f(a,b,x) = atex/D 0<a,0<x. 
(a,b,x) Taba x 


The right tail probability calculated is: 


N 
P(x >a) = f? flab,x) dx. L 


If a/b is 100 or less, two different continued fractions are used to evaluate the significance 
level (References 1 and 3). If a/b is greater than 100, a normal approximation is used 
(Reference 2). 


If a is an even integer and b = 2, the gamma distribution reduces to the chi-square distribution. 
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GAMMA DISTRIBUTION 


REFERENCES 


1. Khovanskii, A. N., The Applications of Continued Fractions and Their Generalizations 
to Problems in Approximation Theory, translated by P. Wynn (1963), Groningen: P. 
Noordhoff. 


2. Peizer, D. B., Pratt, J. W., "A Normal Approximation for Binomial, F, Beta, and Other 
Common Related Tail Probabilities, |”, Journal of the American Statistical Association, 
63 (1968), pp 1416—1456. 


3. Wall, H. S., Analytic Theory of Continued Fractions, New York: D. Van Nostrand, 1948. 


VARIABLES 


VARIABLE USAGE 
PN A Parameter a 
| A1 Scratch 
A2 Scratch 
n B Gamma/2 
B1 Scratch 
n B2 Scratch 
t C1 Scratch 
C2 Parameter b 
1 F Value of gamma 
l | Counter 
N N Gamma 
P1 Scratch 
T Scratch 
1 X Scratch 
1 
"^ 
n Plot 50: Stat Vol 1 l 1 -1 1 


DISTRIBUTIONS 


C 
C C GCG G E: 


GAMMA DISTRIBUTION 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 4; press RETURN. 
4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 5; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the 1/0 light goes out, execute the program. 


C GG Z ES 


C 


CC 


PROGRAM EXECUTION 


1. When the message A, B, GAMMA is displayed, enter the values and press RETURN. || 
Separate value entries with commas. 


2. The values just entered and the significance level are displayed. i 


3. Other values may be entered for further calculations, or the BREAK key may be pressed 
twice to exit from the program. 


C 


EXAMPLE d 


GAMMA DISTRIBUTION 


A B GANNA 58,2,62.3 i 
Az58 B22 PCGAMMA> 62.3 2=0.11372 


A,B, GAMMA 4.5,1,2.3 
824.5 Bai PCGAMMAD 2.3 )=@,86769 


AB GAMMA 12.3,2.6.1.9 | 
0212.1 B=2.6 PCGAMMA> 1.9 )=0.99836 (O 
As B, GAMMA 

1 -1 2 Plot 50: Stat Vol 1 
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NORMAL DISTRIBUTION 


DESCRIPTION 
This program calculates one-tail and two-tail significance levels for a normal distribution. 


The x-value is the only required input. One-tail and two-tail significance levels are displayed 
with the x-value. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


A normal approximation is used with a maximum error of 7.5 x 10°. Increased accuracy is 
available with the High Accuracy Normal program described in this section. 


METHODS 


The approximation used by this program is equation 26.2.17 of the reference. 


REFERENCE 


Abramowitz, M. and Stegun, I. A., editors, Handbook of Mathematical Functions with 
Formulas, Graphs, and Mathematical Tables, Washington, D. C: National Bureau of Standards, 
Applied Mathematics Series No. 55, Government Printing Office. 


VARIABLES 


VARIABLE USAGE 


Scratch 
Scratch 
Normal deviate 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 
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NORMAL DISTRIBUTION 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 5; press RETURN. 
When the I/O light goes out, execute the program. 


> 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 6; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the |/O light goes out, execute the program. 


C 


PROGRAM EXECUTION 
1. The displayed message is the letter X. Enter the X-value and press RETURN. 


Co Ch 


2. The X-value and the significance level are displayed. 


3. Additional values may be entered, or the BREAK key may be pressed twice to 
terminate the program. 


C. C 


EXAMPLE 
NORMAL DISTRIBUTION (APPROXIMATE) U 
X 1.2 | 
X=1.2 ONE TAIL=0.11507 TWO TAIL=0. 23814 
X -2.1 
Xz-2.1 ONE TAIL=0.017864 TWO TAIL=0, 035728 UJ 
X 1.98 
X=1.98 ONE TAIL=0.023852 TWO TAIL=0.047704 VÀ 
) U 


1 = 1 4 Plot 50: Stat Vol 1 U 
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BETA DISTRIBUTION 


DESCRIPTION 


This program calculates right tail probability for the beta distribution with parameters a and 
b. For right tail probabilities between .0001 and .9999, accuracy of 5 digits or more is pro- 
vided when both parameters are 500 or less. When both parameters are greater then 500, 3 
or more digit accuracy is calculated. The a and b parameters and the X-value required as in- 
put are returned with the right tail probability as output. 


HARDWARE REQUIREMENTS 


- A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


Both parameters should be 0.1 or greater. 


METHODS 
The density function of the beta distribution is: 


a—1 (1-x)b—1 
g(a,b,x) aa a>0,b>0,and0<x<1 
B(a,b) 


where B(a,b) is the complete beta function. 


The right tail probability for parameters a and b and value X is: 


M g(a,b,t) dt. 


When either a or b are 500 or less, a continuous fraction is used to evaluate right tail pro- 
bability (Reference 1 and 2). When both are greater than 500, a normal approximation is 
used (Reference 3). 


REFERENCES 


1. Aroian, T. A., "Continued Fractions for the Incomplete Beta Function", Annals of 


Mathematical Statistics, 12 (1941), pp 218—223. 
2. Kopitzke, R. W., “Unpublished Notes”. 


3. Peizer, D. B., Pratt, J. W., "A Normal Approximation for Binomial, F, Beta, and Other 
Common, Related Tail Probabilities, I", Journal of the American Statistical Association, 
63, (1968), pp 1416—1456. 
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VARIABLES 


VARIABLE USAGE 

A Parameter A L 
Al Scratch 

A2 Scratch 
B Parameter B i 
B1 Scratch 

B2 Scratch 

| Scratch U 
Q Value of beta 

T Scratch 

U Scratch L 
V1 Scratch 

X Scratch 

Y4 Scratch U 


C 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


ES uL. E 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 6; press RETURN. 
4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 7; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the I/O light goes out, execute the program. 
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BETA DISTRIBUTION 


PROGRAM EXECUTION 
. A, B, X is displayed. Enter the values separated by commas, and press RETURN. 
. The input values and the calculated right-tail probability are displayed. 


. Other values may be entered in the same manner, or the program may be terminated 
by pressing the BREAK key twice. 


BETA DISTRIBUTION 


AsBsX 50,58,.5 
A=30 B=50 PC(X/0.5)2=0,5 


As Bi X 24,26.35,.8 
fin24 B=26.35  P(450.8)22,.0606E-? 


&,B,X 1.58,4.36,.12 
fx1.58 — Bz4.36 —P(4500.12220.78147 


LITT 
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POINT BINOMIAL DISTRIBUTION 


DESCRIPTION 


This program calculates the individual terms of the binomial distribution for input para- 
meters N, R, and P. N denotes the number of trials; R, the number of successes; and P, 
the probability of success. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 
P(X = R) must be greater than 10 °°°., 


METHODS 


For each set of parameters, the quantity calculated is 


NI 
P(X =R) = ————— PR(1_p)N-R 
(N=R) ! R! 


C 
C C Ef. L- 


If N is less than 50, and R is less than 20, exact factorials are calculated. Otherwise, 
Sterling's approximation is used. 


C 


Probabilities are calculated with a relative error less than 10 *. 


REFERENCE 
Dixon, W. J., Massey, F. J., /ntroduction to Statistical Analysis, New York: McGraw-Hill, 
1957. 
U 
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VARIABLES 


VARIABLE USAGE 


A1 Scratch 

B1 Scratch 

N Number of trials 

P Probability of success 
Number of successes 

X1 Scratch 

X2 Scratch 


PROGRAM LOADING 


us | 


FY This program may be loaded automatically from the directory, or manually from the 
keyboard. 


Lea; 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 7; press RETURN. 


md od 


4. When the I/O light goes out, execute the program. 


— 
e. 


od 22 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 8; press RETURN. 
3. Enter OLD: press RETURN. 
4 


st 


. When the I/O light goes out, execute the program. 
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POINT BINOMIAL DISTRIBUTION 


PROGRAM EXECUTION 


NUMBER OF TRIALS, NUMBER OF SUCCESSES, PROB OF SUCCESS is displayed. 
Enter the values separated by commas, and press RETURN. 


T 


EXAMPLE 


1-20 


. Input values and calculated probability are displayed. 


. Repeat the program execution process for other values, or press BREAK twice to exit 
from the program. 


C 


E op. Ee ES: d Lo Lo Cs 


CC 


POINT BINOMIAL DISTRIBUTION 


NUMBER OF TRIALS,NUMBER OF SUCCESSESSPROB OF SUCCESS 56,23,.2 
N=56 P20.2  P(Xz23021.684E-4 


NUMBER OF TRIALS, NUMBER OF SUCCESSES,PROB OF SUCCESS 10,9,.8 
N=10 P20.8  P(X255»0.026424 


NUMBER OF TRIALS, NUMBER OF SUCCESSES,PROB OF SUCCESS 24,23,.2 
N=24  PxQ.2  P(X-230:1.6186E-15 


HUMBER OF TRIALS,NUMBER OF SUCCESSES.PROB OF SUCCESS 


C 
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SINGLE TERM DISTRIBUTION 


DESCRIPTION 


For the binomial distribution, this program calculates the probabilities for each term be- 
tween designated lower and upper limits. The accumulative and right tail probabilities are 
also calculated. 


Required input values are the number of trials, lower and upper term limits, and the pro- 
bability of success. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


For probability of success P, and number of trials N, (1—P)N must be greater than 10°°°. 


The number of trials should be 50 or less. For more than 50 trials, the Accumulative Bi- 
nomial Distribution program should be used. 


METHODS 


If X has a binomial distribution with parameters shown above, the probabilities calculated 
are as follows: 


N 
single term; P(X =r) = pr (1-P)N-", 
r 
j 
accumulative; P(X <j) =  P(X=i) 
i=0 


right tail; P(X >j) = 1—P(X <j). 


P(X = 0) is first calculated. Probabilities for X greater than O are then calculated using the 
recursion relation: 


_ (N—j+1) P 
j 1-P 


P(X = j) P(X = j—1) for j 1,2,...,N. 
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REFERENCE 


Dixon, W. J., Massey, F. J., Introduction To Statistical Analysis, New York: McGraw-Hill, 
1957. 


VARIABLES 


VARIABLE USAGE 


| Scratch 

N Number of trials 

P Probability of success 

P1 Scratch 

P2 Scratch q 
P3 Scratch 

P4 Scratch e 
R Number of successes V J 


PROGRAM LOADING 
This program may be loaded automatically from the directory, or manually from the 
keyboard. 


EE 


p 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory wil! be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 8; press RETURN. 
4. When the I/O light goes out, execute the program. 


U 
MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. L 
2. Enter FIND 9; press RETURN. 
3. Enter OLD; press RETURN. V ; 
4. When the I/O light goes out, execute the program. H 
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N 
25 
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PROGRAM EXECUTION 


1. NUMBER OF TRIALS, LOWER LIM, UPPER LIM, PROB OF SUCCESS appears on 
the display screen. Enter values for the number of trials, lower limits, upper limits, and 
probability of success. A comma must separate the values. Press RETURN. 


DISTRIBUTIONS 


SINGLE TERM DISTRIBUTION 


2. Input values and calculated probabilities are displayed. 


INGLE TERM BINOMIAL 


POGY) 
8.8312500 
0. 1362300 
8. 3125000 
9. 3123000 
8. 1562500 
0. 8312588 


POKRY) 
0.0312500 
0. 1875008 
8. 5000000 
8.8125980 
0. 9697588 
1. 0000090 


S 
N=5 Pæ0.3 PROBABILITIES CALCULATED FROM @ TO 5 SUCCESSES 


P(X>Y) 
0. 9687500 
8.812580 
9. 9000000 
8. 1875000 
8. 8312508 
0. 0000060 


SINGLE TERM BINOMIAL 


PSY) 
9.8916360 


POD YO 


09.0979000 


3. To re-execute this program, type RUN, press RETURN, and enter new data values. 


Find the probability of success for all possible values for a binomial distribution with 5 
trials and .5 probability of success. 


For a binomial distribution with probability of success .3, find the probability of 10 suc- 
cesses in 25 trials. 


N=25 P=0.3 PROBABILITIES CALCULATED FROM 10 TO 10 SUCCESSES 


POP) 
9.9022008 
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SINGLE TERM DISTRIBUTION 


For a binomial distribution with probability of success .2, find the probability of each suc- 


cess between 90 and 110 for 500 trials. 


N=500  P=0.2 PROBABILITIES CALCULATED FROM 98 TO 118 SUCCESSES 


1-24 


SINGLE TERM BINOMIAL 


P(X#Y) 
0.0244260 
8, 0275129 
8. 0305781 
8. 8335373 
8. 0363024 
0. 0387862 
0. 0409073 


0. 0321654 
à. 8292616 
8. 0263086 
8.8233788 


PCX<2Y) 


8. 7329694 
0.7679105 
8. 8000759 
09.8293375 
6, 8556461 
8,8798249 


PCX>Y) 
8.8562972 
8,8297844 
0. 7982062 
8,7646689 
8,7283666 
9,6895803 
8,6486739 
8. 6960788 
6,5622893 
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C 


[= Es 


C 


C 


C Lã 


Plot 50: Stat Vol 1 


2 


ue. Id 


ee Nesen a 3 .—3 SI 
2 


DISTRIBUTIONS 


ACCUMULATIVE BINOMIAL DISTRIBUTION 


DESCRIPTION 


This program calculates the accumulative binomial distribution only. The input parameters 
are the number of trials and successes, and the probability of success. Output is the accu- 
mulative probability of successes. 


When the number of trials is 500 or less, a continued fraction is used to evaluate the pro- 
bability (References 1 and 2). Minimum accuracy is five digits. If trials minus successes plus 
1 is greater than 500, a normal approximation is used and minimum accuracy is two digits 
(Reference 3). 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


This program is intended primarily for use when the number of trials exceeds 50. The single 
term binomial distribution program should be used when the number of trials is 50 or less. 


METHODS 


When N represents the number of trials, R the number of successes, and P the probability 
of success, the binomial distribution is defined by: 


N! 
P(X=R)= ——— — PÅ (1-pjN- R. 
(N-R)!R! 
The quantity calculated is: 


2 P(X =i) - 


For the parameters above, P(X > R) is a beta with parameters a = R and b = N-R+1, 
evaluated at X = P. The methods used are the same as in the beta program. 
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REFERENCES 


1. Aroian, T. A., "Continued Fractions for the Incomplete Beta Function", Annals of 
Mathematical Statistics, 12 (1941), pp. 218—223. [| 


2. Kopitzke, R. W., “Unpublished Notes.” 


3. Peizer, D. B., Pratt, J. W., “A Normal Approximation for Binomial, F, Beta, and Other 
Common Related Tail Probabilities, 1’ Journal of the American Statistical Association, 
63 (1968), pp. 1457—1483. 
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VARIABLES | 
VARIABLE USAGE U 
A R+1 
A1 Scratch å 
A2 Scratch vw i 
B N—R 
B1 Scratch L 
B2 Scratch 
| Counter 
N Number of successes E 
P Probability of success 
Q Value of beta lU 
R Number of trials 
T Scratch : 
U Scratch i 
V1 Scratch 
X Scratch 
YA Scratch U 
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ACCUMULATIVE BINOMIAL DISTRIBUTION 
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PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 9; press RETURN. 
4. When the !/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 10; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the I/O light goes out, execute the program. 


5 


PROGRAM EXECUTION 


1. NUMBER OF TRIALS, NUMBER OF SUCCESSES, PROB OF SUCCESS is displayed. 
Enter the respective values separated by commas. Press RETURN. 


2. The accumulative probability will be displayed with the input values. 


3. Other values may be entered to execute the program again, or an exit may be called 
by pressing the BREAK key twice. 


EXAMPLE 


ACCUMULATIVE BINOMIAL DISTRIBUTION 


NUMBER OF TRIALS,NUMBER OF SUCCESSES,PROB OF SUCCESS 590,97,.2 
N=308 P=0.2 PCX<=97)=8, 39392 


NUMBER OF TRIALS, NUMBER OF SUCCESSES,PROB OF SUCCESS 100091, 500098, .5 
N=100001 P=0.5 P(X<250000)20.5 


NUMBER OF TRIALS,NUMBER OF SUCCESSES,PROB OF SUCCESS 50000, 2580,.9 
N=50008 P=0.9 P(X(22900520 


NUMBER OF TRIALS, NUMBER OF SUCCESSES, PROB OF SUCCESS 
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DESCRIPTION 


This program calculates the probability that a hypergeometric random variable has the value 
Y. The random variable represents the number of defectives in sample size M of population 
size N. K represents defectives in the population, and Y represents defectives in the sample. 


The values of N, M, K, and Y are required input. 


All probabilities are calculated with minimum six digit accuracy. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 
P(X = Y) must be greater than 103009. 


METHODS 


All factorials are calculated using Stirlings approximation. 


The quantity calculated is: 


hee 


in 


P(X =Y)= 


REFERENCE 


Abramowitz, M., and Stegun, I. A., editors, Handbook of Mathematical Functions with 
Formulas, Graphs, and Mathematical Tables, Washington, D. C: National Bureau of Standards, 


Applied Math Series No. 55, Government Printing Office. 
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VARIABLES 
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PROGRAM LOADING 


VARIABLE USAGE 


n 
1 
A1 Scratch 
n B1 Scratch 
K Number of defectives in the 
population 
f Sample size 
| N Population size 
q X Number of defectives in the 
Å sample 
X1 Scratch 
q X2 Scratch 
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POINT HYPERGEOMETRIC DISTRIBUTION 


This program may be loaded automatically from the directory, or manually from the 


keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 


1. Insert the program tape. 


THE PROGRAM NUMBER YOU WANT. 
. Enter 10; press RETURN. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 

. Enter FIND 11; press RETURN. 

. Enter OLD; press RETURN. 


RR WN 
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4. When the I/O light goes out, execute the program. 


. When the 1/0 light goes out, execute the program. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
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PROGRAM EXECUTION 


1. 


LOT SIZE is displayed. Enter the population value and press RETURN. 


2. SAMPLE SIZE is displayed. Enter the value and press RETURN. 


Co 


Oo oc 


. NUMBER OF DEFECTIVES IN THE LOT appears on the screen. Enter the value; 


press RETURN. 


. NUMBER OF DEFECTIVES IN THE SAMPLE is displayed. When the value is entered, 


press RETURN. 


. The probability is displayed with the values which were entered. 


. Steps 1 through 5 may be repeated for other values, or the program may be terminated 


by pressing the BREAK key twice. 


EXAMPLE 


POINT HYPERGEOMETRIC DISTRIBUTION 


LOT SIZE 58 
SAMPLE SIZE 24 
NUMBER OF DEFECTIVES IN THE LOT 28 
NUMBER OF DEFECTIVES IN THE SAMPLE 16 
N= 50 Ms 24 K= 20 P(X* 16)= 0.221045 


LOT SIZE 560 

SAMPLE SIZE 108 

NUMBER OF DEFECTIVES IN THE LOT 25 

NUMBER OF DEFECTIVES IN THE SAMPLE 12 
N= 560 M= 106 K= 25 P= 13507 


ce 
ort] 
ce 
[Xx] 
@ 
o 
ma 


LOT SIZE 
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SINGLE TERM HYPERGEOMETRIC DISTRIBUTION 


DESCRIPTION 


This program calculates probabilities for each value of random variable X between the upper 
and lower limits of Y. X represents the number of defectives in the sample and follows the 
hypergeometric distribution. Y limits are determined at execution time. 


Probabilities are calculated with 10 digit minimum accuracy. 


Population size (N), sample size (M), the number of defectives in the lot (K), and the upper 
and lower limits to the number of defectives in the sample (Y), are the required input values. 


Output probabilities are expressed as X = Y, X <= Y, and X > Y. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


When accumulative distribution (X <= Y) is the only output desired, the Accumulative 
Hypergeometric Distribution program should be used. 


METHODS 


The hypergeometric function is defined as: 


mr) lg) 


P(X =R) =| M-R RI. 


(Ml 
M 
Let R = maximum (0,M+K—N). P(X = R) is first calculated and the recursion relation 


(K—J+1) (M—J+1) 
P(X=J) = (N-K-MH) ^ J *P(X=J-1) 


is used to calculate other probabilities. 
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REFERENCE 


Mood, A., Graybill, F., Boes, D., /ntroduction to the Theory of Statistics, New York: 
McGraw-Hill, 1974. 


ES fo ta 


VARIABLES 


USAGE 


Number of defectives in the population 


VARIABLE 


Sample size 

Population size 

Scratch 

Number of defectives in the sample 
Scratch 

Scratch - | | 
Scratch t 
Scratch 


Scratch 7 


Scratch 


E EB © fe 


[-— E 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 
= 1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 11; press RETURN. 
4. When the I/O light goes out, execute the program. 


U— 


C 


Plot 50: Stat Vol 1 


CZ 


Z4 


2 


—J 


) 


) 


— eo X. LJ 


n | 


DISTRIBUTIONS 


SINGLE TERM HYPERGEOMETRIC DISTRIBUTION 


MANUAL LOADING FROM THE KEYBOARD 


PWN 


. Insert the program tape. 

. Enter FIND 12; press RETURN. 

. Enter OLD; press RETURN. 

. When the I/O light goes out, execute the program. 


PROGRAM EXECUTION 


1. 


LOT SIZE appears on the display screen. Enter the lot (population) size and press 
RETURN. 


2. SAMPLE SIZE is now displayed. Enter the sample size; press RETURN. 
3. When NUMBER OF DEFECTIVES IN THE LOT is displayed, enter the value and 


press RETURN. 


. CALCULATE PROBABILITIES FOR NUMBER OF DEFECTIVES IN SAMPLE 


> = is displayed. This value is the lower limit for probability calculation. Enter the 
value and press the RETURN key. 


. AND FOR NUMBER OF DEFECTIVES IN SAMPLE < = appears. Enter the upper 


probability limit and press RETURN. 


. Probabilities for the number of defectives in the sample which are within the limits 


are printed. 


. To repeat execution of this program, type RUN, press RETURN, and enter new data 


values. 
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SINGLE TERM HYPERGEOMETRIC DISTRIBUTION 


EXAMPLE 1 


For population 40, sample 23, and 12 defectives in the population, calculate probabilities 


for all possible defectives in the sample. 


Set the lower limit at zero and the upper limit high enough to calculate al! values. In this 
example, 12 is the highest possible upper limit value. After all input values have been entered, E 


calculations will be printed when the RETURN key is pressed. 


SINGLE TERM HYPERGEOMETRIC DISTRIBUTION 
FOR LOT SI2E= 40 SAMPLE SIZE= 23 NUMBER OF DEFECTIVES IN LOT = 12 
CALCULATED Au MER DEFECTIVES IN SAMPLE BETWEEN @ AND 12 


Y PCE2Y) PCX<2¥) PODYO 
42 23 0 6.000081 8.000001 8.999999 
48 23 i 0.006051 20.000052 8.999948 
40 23 2 80.000881 6.800932 6.999067 
40 23 3 8.807786 4.009639 8.991361 
42 23 4 90.838521 8.047178 8.952838 
40 23 5 0.117133 60.164303 — 0.83569? 
40 23 6 0,223619 6.387921 6.612879 
40 23 ? 8.271526 8.659456 — 8.340544 
40 23 8 0.200874 0.968328 — 0.121670 
48 23 9 0.039464 4.967794 9.032266 
40 23 18 8.027850 9.993644 0.884356 
40 23 11  e.8804114 6.929758 — 0,500242 
48 23 12 0.008242 — 1.008000 0.800800 
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SINGLE TERM HYPERGEOMETRIC DISTRIBUTION 


EXAMPLE 2 


The lot size is 500, the sample size is 50, and there are 100 defectives in the lot. Calculate 
probabilities for all possible number of defectives in the sample with non-zero probabilities. 


Since the non-zero probability limits are unknown, arbitrary limit assignments are made. In 
“ this example, limits of 10 and 40 were chosen. When the RETURN key is pressed, the calcu- 
lated results will appear on the screen. 


SINGLE TERM HYPERGEOMETRIC DISTRIBUTION 
FOR LOT SIZE= 500 SAMPLE SIZE= 50 NUMBER OF DEFECTIVES IN LOT = 108 
CALCULATED FOR NUMBER DEFECTIVES IN SAMPLE BETWEEN 19 AND 48 
N M Y POiSY) POS Y) POUDY) 
0.147368 6.585148 0.414852 


3 ud. 3 
g 


500 58 11 0.133600 20.718748 0.281252 

$00 5e 12 9.106751 80.825499 0.174501 

580 50 13 0.075646 80.901145 6.698855 

$00 38 14 0.047784 40.543928 9.051872 

388 38 15 8.827021 5.375943 6.024051 

380 56 16 0.01372? 9.989676 X 0.019224 

580 50 1?  0.0%6234 0,335960 0. 004040 

es 500 50 19 60.082595 6,9985589 8.001442 
j| FON 588 38 19 0.000972 3.999531 0.000459 
509 58 26 80.000330 4.999861 . 0.800139 
500 56 21  $8.900102 0,999962 4,069438 

588 58 22  UU.000020 «= 99991 — 0,008005 

n 500 50 ei 9.600007 .999398 — 8.000802 
| 380 38 e 0.020862 agagap — 3.900085 
S80 38 25 à, 846888 Hung X 3.0000n0 


289000 X G.000C80 
002090 X 6.000008 
000006 X 0.000008 
860080 6.000086 
090000 6.008660 
568 50 32 80.000008 1.000800 9.080000 
$09 50 33 0.000000 1.086009 0.080080 
300 S8 34 0.000008 1.000000 8.680088 
500 38 35 80.000060 1.088080 0.080000 
568 56 36 80.080000 1.000000 0.000080 
588 50 7 0.000000 1.006800 8.080088 
500 38 38 0.000000 1.908008 6.000000 
$80 58 39 0,000000 1.090000 6.800880 
500 sø 40 6.000008 1.000000 6.880000 


"900900 0.880808 


«a 
e 
© 
cn 
© 
e 
Lnd 
e 
e 
e 
Q 
oa 
& 
«e 


The table printed above indicates that the lower limit (10) is not low enough and that an 
upper limit of 25 is correct. 
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When limits are set between 0 and 25, non-zero probabilities are shown as: 


SINGLE TERM HYPERGEOMETRIC DISTRIBUTION 


FOR LOT SIZE= $88 SAMPLE SIZE= 50 NUMBER OF DEFECTIVES IN LOT = 108 


CALCULATED å BARER DEFECTIVES IN SAMPLE BETWEEN Å np 25 


Y 


OG OD 74 € CR å CJ 03 à OO BIH HS Gå ND e CD 


Put eu br Dodo puto pode nt pt put qua 


28 


8. 008090 


PCX<=Y) 
8.088007 
0.080112 
0. 000835 
0. 004044 
0.01437? 
B. 040034 
0.691534 
B. 176687 
9.292587 
0.437781 
8.585148 
» 718748 
. 825499 
. 901145 
. 948928 
2975949 
«989676 
. 995360 
. 998558 
. 993531 
. 939861 
999962 
6.999991 
8.999998 
1. 868008 
1. 000090 


OD € On G2 (S Orage 


PX 
ð. 599993 
0.999888 
0.999165 
8.995956 
8, 985623 
0.959916 
8.908466 
8.823313 
0.704413 
9.562219 
8.414852 
0.281252 
0.174581 
0.099955 
8.451872 
9.024851 
9.810224 
9.004848 
8.001442 
8.608463 
8.600139 
a. 000038 
9. 008605 
9. 000002 
8.000000 
8.608000 
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ACCUMULATIVE HYPERGEOMETRIC DISTRIBUTION 


DESCRIPTION 


This program calculates the probability that Y is less than or equal to R. Y is a random vari- 
able following the hypergeometric distribution, and R is a selected number of defectives in 
the sample. 


The accumulative probability is calculated with a minimum of 10 significant figures of 
accuracy. 


Lot size (N), sample size (M), number of defectives in the lot (K), and number of defectives 
for which probability will be calculated (R), are required input. 


HARDWARE REQUIREMENTS 


A Tektronix 4050-Series Graphic System is required. 


METHODS 


Sterlings approximation is used to calculate P(Y = R). 


The recursion relation is 
(Y—j—1) (N-M-K+Y—j—1) 


P(Y = j—1) = Se ee | Se «P(Y=j) 
(M-Y-1—j) (K-X*1-j) | 


is used to calculate lower order probabilities. 


REFERENCE 


Mood, A., Graybill, F., Boes, D., /ntroduction of the Theory of Statistics, New York: 
McGraw-Hill, 1974. 
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ACCUMULATIVE HYPERGEOMETRIC DISTRIBUTION 


VARIABLES 


Co C 


VARIABLE USAGE 
Scratch 


E 


Al Scratch 

A2 Scratch 

A3 Scratch | 
B1 Scratch 
B2 Scratch || 
K Number of defectives in population 

M Sample size i 
N Population size 

X Number of defectives in sample | 
X1 Scratch u 
X2 Scratch 


q 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


[T 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 12; press RETURN. 
4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
. Enter FIND 13; press RETURN. 
. Enter OLD; press RETURN. 
. When the I/O light goes out, execute the program. 
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DISTRIBUTIONS 


n ACCUMULATIVE HYPERGEOMETRIC DISTRIBUTION 


4 PROGRAM EXECUTION 
! 1. When LOT SIZE is displayed, enter the population value and press RETURN. 


2. SAMPLE SIZE indicates that the sample value is needed. Enter the value and press 
RETURN. 


3. NUMBER OF DEFECTIVES IN LOT appears. Type the value and press RETURN. 


4. CALCULATE PROBABILITY OF NUMBER OF DEFECTIVES IN SAMPLE <= 
requests the final required value. Enter the value; press RETURN. 


c 


5. Input values and calculated probability will appear on the screen. 


6. Repeat the program execution sequence for other values, or press BREAK twice to 
exit from the program. 


EXAMPLE 


For lot size 50, sample size 23, and 12 defectives in the lot, find the probability that the 
PN sample will contain 6 defectives or less. 


cJ es ese, 


23 


ACCUMULATIVE HYPERGEOMETRIC DISTRIBUTION 
LOT SIZE 50 
SAMPLE SIZE 23 
NUMBER OF DEFECTIVES IN LOT 12 
VENNE Tel) OF NIMES "Mr E IN SAMPLE <= 6 


3 


Ne 23 K ¿a 6)* 8,742791 
1 LOT SIZE 
7 
i 
» 
I 
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1-40 


This program calculates probabilities for a random variable following the Poisson distribution 
with rate parameter (mean) lambda. Probabilities are calculated for each value between 
specified lower and upper limits. Accumulative and right tail probability are also calculated. 


U 
DESCRIPTION | | 
i 


The values of the mean and the lower and upper limits must be entered. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. | 


PROGRAM LIMITS i 


The rate parameter must be positive. 


METHODS 
With non-negative integer r and rate parameter m, the Poisson distribution is defined by KJ Ll 
P(X =r) = ec m" T 
r! | 
REFERENCE 


Mood, A., Graybill, F., Boes, D., /ntroduction to the Theory of Statistics, New York: 
McGraw-Hill, 1974. 


E. ES 


VARIABLES 


E 


VARIABLE USAGE 


Accumulative probability 
Exp (— lambda) 


CGC ES. 


Counter 
Lower limit 
Upper limit 
Scratch 
Scratch 


C 
C fI OC 


Mean lambda 
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PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


ot 23 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 13; press RETURN. 
4. When the I/O light goes out, execute the program. 


23 53 


MANUAL LOADING FROM THE KEYBOARD 

. Insert the program tape. 

Enter FIND 14; press RETURN. 

Enter OLD; press RETURN. 

When the 1/0 light goes out, execute the program. 


) 


>ON 


a LL 3 


z3 


PROGRAM EXECUTION 
1. LAMBDA is displayed. Enter the value of the mean and press RETURN. 
2. CALCULATE PROBABILITIES: 


FROM LOWER LIMIT appears on the screen. Enter the smallest value for which 
probability is to be calculated. Press RETURN. 


. TO UPPER LIMIT is displayed. Enter the largest value for probability calculation and 
press RETURN. 


4. Probabilities are displayed for all values within the limits. 


5. Program execution may be repeated by typing RUN, pressing RETURN, and entering 
other data values. 


QD 


2 
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SINGLE TERM POISSON DISTRIBUTION 
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EXAMPLE 1 


For a Poisson variable with mean 4, calculate probabilities for values from zero to 10. 


DW 00 ^40 UT & GIRO v 


vat 


EXAMPLE 2 


8.005292 


SINGLE TERM POISSON DISTRIBUTION 


PCX>Y) 
8. 


PCX<2Y) 
0.018316 
0.091578 
0.238103 
8.433478 
0.628837? 
8.785136 
0.889326 
8.948366 
6.978£37 
0.931868 
8,997168 


981684 


For a Poisson variable with mean 6.2, calculate the probabilities of values between 2 and 18. 


LANBDA = 6.2 


Y 


A» U AD bet 0930 OO ^40 CA > OG PI 


pt forde Quo Jonte 


tua me 
Hu 
Sa 
sa O 
DA 
Sa 


te 
e 
d 
co 
wo 
^J 


0.000058 


SINGLE TERM POISSON DISTRIBUTION 
PODY) 


POXCS YO 


«339315 
«995973 
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DISTRIBUTIONS 


HIGH ACCURACY NORMAL DISTRIBUTION 


DESCRIPTION 


This program calculates right tail probability for the normal distribution. Probabilities are 
computed with at least 10 digit accuracy. A normal deviate value is the only required input. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


Normal deviate values must be positive. 


For normal deviates greater than 20, the log to base 10 of the normal probability is printed. 


METHODS 


An infinite series is used to calculate probability for normal deviates of 3.5 or less. When 
the deviate is greater than 3.5, Gauss's continued fraction is used. 


REFERENCE 


Abramowitz, M., and Stegun, |. A., editors, Handbook of Mathematical Functions with 
Formulas, Graphs, and Mathematical Tables, Washington, D. C: National Bureau of Standards, 
Applied Math Series No. 55, Government Printing Office. 


VARIABLES 


VARIABLE USAGE 


A1 Scratch 
A2 Scratch 
B1 Scratch 
B2 Scratch 
| Counter 
X Absolute value of normal deviate 
Y Scratch 
Y1 Scratch 
Z l Scratch 
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HIGH ACCURACY NORMAL DISTRIBUTION V I 
PROGRAM LOADING | i 
This program may be loaded automatically from the directory, or manually from the L 

keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. | 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 14; press RETURN. 
4. When the I/O light goes out, execute the program. 


i 


(On 


MANUAL LOADING FROM THE KEYBOARD 


1. Insert the program tape. 1 
2. Enter FIND 15; press RETURN. 

f 
3. Enter OLD; press RETURN. \ / i 


4. When the I/O light goes out, execute the program. 


PROGRAM EXECUTION 


1. X is displayed on the screen. Enter the value for the positive normal deviate; press [| 
RETURN. 

2. The right tail probability is displayed. i 

3. Another value may be entered, or the program may be terminated by pressing the i 


BREAK key twice. 
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HIGH ACCURACY NORMAL DISTRIBUTION 


EXAMPLE 


HIGH ACCURACY NORMAL DISTRIBUTION 
X 


X 1.2 

PCX>1.2) =0,115069670222 
X 2.6 

P(X>2.6) 20.0046611880236 


X 9.4 
PODS9,4) = 2,728153571E-21 


X -1 
P(X>1) 20.158655253932 


P oll JJ.  —J. l3 3  —3J 


X 27 
LOG(P(X>27)) = -160.131286062 
A 


23 Z3 
D 


ba ss Gene n À 


23 


) 
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TABLED VALUES 


F-TABLED 


DESCRIPTION 


This program calculates the tabled value for an F-Distribution with degrees of freedom N 
and D, and right tail probability P. The program is the inverse of the F-Distribution pro- 
gram in the DISTRIBUTIONS section. 


Accuracy is calculated to a minimum of 6 digits for significance levels between .0001 and 
.9999. 


Degrees of freedom and significance level (right tail probability) are required input. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


Degrees of freedom must be greater than .1, and the significance level must be between 
.0001 and .9999. 


METHODS 


The F-distribution is transformed to a beta distribution and a tabled beta estimate is calcu- 
fated. A quasi-linear interpolation scheme uses this estimate to solve for the tabled beta 
value. The tabled beta value is then transformed to a tabled F value (Reference 2). 


REFERENCES 


1. Aroian, T. A., "Continued Fractions for the Incomplete Beta Function”, Annals of 
Mathematical Statistics, 12 (1941), pp. 218—223. 


2. Kopitzke, R. W., "Unpublished Notes.'' 


3. Peizer, D. B., Pratt, J. W., "A Normal Approximation for Binomial, F, Beta, and Other 
Common Related Tail Probabilities, |” Journal of the American Statistical Association, 
63 (1968), pp. 1457—1483. 
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TABLED VALUES 


F-TABLED 


VARIABLES 


2-2 


VARIABLE 


USAGE 
Numerator degrees of freedom/2 
Scratch 


Scratch 

Denominator degrees of freedom/2 
Scratch 

Scratch 

Numerator degrees of freedom 
Denominator degrees of freedom 
Counter 

Significance level 

Scratch 

Scratch 

Value of beta function 
Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Approximate root 
Approximate root 
Approximate root 

Scratch 

Scratch 

Scratch 

Scratch 
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TABLED VALUES 


F-TABLED 


—— ^ — 01 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


Z3 


AUTOMATICALLY LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 15, press RETURN. 
4. When the 1/0 light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
. Enter FIND 16; press RETURN. 
. Enter OLD; press RETURN. 
When the 1/0 light goes out, execute the program. 


— a om vn (eee Ves | 


RON 


PROGRAM EXECUTION 


1. NUM DF, DEN DF, SIG LEVEL is displayed. Enter the values and press RETURN. 
Commas should be used to separate values being entered. 


. The tabled value is displayed. 


co 
N 


3. The program may be repeated by entering new values. To exit from the program, press 
BREAK twice. 


EXAMPLE 


F TABLED 


NUM DF, DEN DF SIG LEVEL 2,2,.05 
NUM DF=2 DEN DF=2 SIG LEVEL=0. 805 TABLED F219 


NUM DF, DEN DF, SIG LEVEL 25,63,.81 
NUM DF=25 DEN DF=63 SIG LEVEL=8.01 TABLED F=2.08228446 


NUM DF, DEN DF SIG LEVEL 253,63, . 3999 
NUM DF=25 DEN DF=63 SIG LEUEL=0.9995 TABLED F=0, 238842268 


NUM DF. DEN DF,SIG LEVEL 5080.5000,.6 
NUM DF=5000 DEN DF=5980 SIG LEVEL=8.6 TABLED F=&, 992859511 
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TABLED VALUES 


t-TABLED 


2-4 


E 


DESCRIPTION 


This program calculates the tabled value for a Student's t distribution and is the inverse of 
the t-Distribution program in Section 1. Accuracy is calculated to a 5 digit minimum. 


E 


Input values are degrees of freedom and significance level for a one-tail test. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


Degrees of freedom must be greater than .1. 


E. ES EG Es: 


Significance level must be between .0005 and .5. 


E 


METHODS 


The t-distribution is transformed to a beta distribution, and an initial estimate of the tabled 
beta value is calculated. A quasi-linear interpolation scheme solves for the tabled beta value 
which is then transformed to t-tabled value. 


C 


REFERENCES 


1. Aroian, T. A., "Continued Fractions for the Incomplete Beta Function"; Annals of 
Mathematical Statistics, 12 (1941), pp. 218—223. 


2. Kopitzke, R. W., "Unpublished Notes." 


p 


3. Peizer, D. B., Pratt, J. W., "A Normal Approximation for Binomial, F, Beta, and Other 
Common Related Tail Probabilities, |” Journal of the American Statistical Association, 
63 (1968), pp. 1457—1483. 
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VARIABLES 


VARIABLE USAGE 


Degrees of freedom/2 
Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Degrees of freedom 


Counter 
Significance level 
Significance level 
Scratch 

Point at which to evaluate beta 
Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 
Approximate root 


) 


Approximate root 
Approximate root 
Scratch 

Scratch 

Log of beta function 


o 
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t-TABLED 
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TABLED VALUES 


t-TABLED 
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PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 16; press RETURN. 
4. When the 1/0 light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 17; press RETURN. 


3. Enter OLD; press RETURN. V 


4. When the I/O light goes out, execute the program. 


PROGRAM EXECUTION 


1. DEGREES OF FREEDOM, SIG LEVEL is displayed. Enter degrees of freedom, a 
comma, and significance level. Press RETURN. 


2. Input and t-tabled values are displayed. 


3. Enter other values to repeat execution of the program, or press BREAK twice to exit. 


EXAMPLE 


T TABLED 


DEGREES OF FREEDOM,SIG LEVEL 16,.81 
DEGREES OF FREEDOM=16 SIG LEVEL=0.01 TABLED T22.58348? 


DEGREES OF FREEDOM,SIG LEVEL 17008,.85 
DEGREES OF FREEDOM=1700 SIG LEVEL=0.05 TABLED T21.64575 


DEGREES OF FREEDOM,SIG LEVEL 12..5 


DEGREES OF FREEDOM,SIG LEVEL 12,.2 
DEGREES OF FREEDOM=i2 SIG LEVEL=0.2 TABLED T=0.87261 


DEGREES OF FREEDOM,SIG LEVEL 


E: 
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TABLED VALUES 


CHI-SQUARE TABLED 


DESCRIPTION 


This program calculates the tabled value for a chi-square distribution and is the inverse of 
the chi-square distribution program described in Section 1. 


Degrees of freedom and the significance level for a one-tail test are required input values. 


The tabled value has a minimum of five digits of accuracy. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


The degrees of freedom must be greater than 0.1 and the significance level must be between 
.0005 and .9995. 


METHODS 


An approximate tabled value is obtained using the Wilson-Hilferty approximation (Refer- 
ence 1, equation 26.4.17). A Newton-Raphson iteration scheme is then employed until the 
tabled value is bracketed above and below. A quasi-linear interpolation scheme uses these 
two bracketing points until the tabled value is found. 


REFERENCES 


1. Abramowitz, M., and Stegun, I. A., editors, Handbook of Mathematical Functions with 
Formulas, Graphs, and Mathematical Tables, National Bureau of Standards Applied 
Mathematics Series No. 55, Washington, D. C: Government Printing Office. 


2. Khovanskii, A. N., The Applications of Continued Fractions and Their Generalizations 
to Problems in Approximation Theory, translated by P. Wynn, P. Noordhoff, 
Groningen: 1963. 


3. Peizer, D. B., Pratt, J. W., “A Normal Approximation for Binomial, F, Beta, and Other 
' Common, Related, Trail Probabilities |”, Journal of the American Statistical Association, 


63, 1968, pp. 1416—1456. 


4. Wall, H. S., Analytic Theory of Continued Fractions, New York; D. Van Nostrand, 
1948. 
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CHI-SQUARE TABLED. 
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VARIABLES 


VARIABLE USAGE 


Degrees of freedom/2 
Al Scratch 
A2 Scratch 
B Final root 
B1 Scratch 
B2 Scratch 
C1 Scratch 
| Counter 
N Degrees of freedom 
P Significance level 
P1 Scratch 
P2 Scratch 
T Scratch V 
X Scratch 
X1 Scratch 
x2 Scratch 
X8 Scratch 
Y1 Scratch 
Y2 Scratch 
Z Scratch 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 


1. 
2. 


3. Enter 17; press RETURN. 


Insert the program tape. 


Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. ' 


C 


4. When the I/O light goes out, execute the program. 
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MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 18; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the I/O light goes out, execute the program. 


tgp. IL. 2) 


PROGRAM EXECUTION 


1. DEG FREEDOM, SIG LEVEL is displayed. Enter the degrees of freedom, a comma, 
and the significance level. Press RETURN. 


2. The Chi-square tabled value is displayed with the input values. 


. Repeat the program execution by entering new values, or exit from the program by 
pressing the BREAK key twice. 


EXAMPLE 


CHI SQUARE TABLED 


DEG FREEDOM,SIG LEVEL 59,.05 
FOR DF=50 SIG LEVEL=0.05 CHI SQUARE TABLED=67,5048 


DEG FREEDOM, SIG LEVEL 4,.9999 
FOR DF=4 SIG LEVEL=0.9999 CHI SGUARE TABLED=9.42841848 


DEG FREEDOM, SIG LEVEL 15800,.2 
FOR DF=15808 SIG LEVEL=0.2 CHI SQUARE TABLED=15145.574 


DEG FREEDOM, SIG LEVEL 


=) 
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CHI-SQUARE TABLED 
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TABLED VALUES 


GAMMA TABLED 


2-10: 


DESCRIPTION 


This program calculates the tabled value for a two parameter gamma distribution with right 
tail probability. It is the inverse of the gamma distribution program and calculates the tabled 
value with a minimum five digit accuracy. 


Parameters and right tail probability are required input values and are displayed with the 


calculated tabled value. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 


The ratio of parameter a to parameter b should be greater than O.1. 


Right tail probability should be between .0005 and .9995. 


C 


METHODS 


An approximate tabled value is obtained using the Wilson-Hilferty approximation (Refer- 
ence 1, equation 26.4.17). A Newton-Raphson iteration scheme is used until the tabled 
value is bracketed above and below. A quasi-linear interpolation scheme uses these two 
bracketing points to find the tabled value. 


REFERENCES 
1. Abramowitz, M., and Stegun, |. A., editors, Handbook of Mathematical Functions 
With Formulas, Graphs, and Mathematical Tables, National Bureau of Standards 
Applied Mathematics Series No. 55, Washington, D. C: Government Printing Office. 


E— 


2. Khovanskii, A. N., The Applications of Continued Fractions and Their Generalizations 
to Problems in Approximation Theory, translated by P. Wynn, 1963, Groningen: 
P. Noordhoff. 


3. Peizer, D. B., Pratt, J. W., “A Normal Approximation for Binomial, F, Beta, and Other 
Common, Related, Tail Probabilities |", Journal of the American Statistical Association, 


63, 1968, pp. 1416—1456. U 


4. Wall, H. S., Analytic Theory of Continued Fractions, New York: D. Van Nostrand, 
1948. me 
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GAMMA TABLED 


VARIABLES 


Mis 2 0523 3 


VARIABLE USAGE 


n Parameter a 
A1 Scratch 

n A2 Scratch 

! B Final root 
B1 Scratch 

A B2 Scratch 
C1 Scratch 


Parameter b 


-3 
Q 
N 


Counter 


Parameter a 


—J 
z 


P Right tail probability 
a PN P1 Scratch 
Fi P2 Scratch 

T Scratch 
1 X Scratch 


Approximate root . 


Approximate root 


23 
S 


X8 Scratch 

n Y1 Scratch 
i Y2 Scratch 
Scratch 


va IR 11 
N 


2 
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GAMMA TABLED 
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PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 18; press RETURN. 
4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 19; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the I/O light goes out, execute the program. 


PROGRAM EXECUTION 


1. A, B, RIGHT TAIL PROB is displayed. Enter the values (separated by commas) and 
press RETURN. 


2. The input values and the tabled value are displayed. 


3. Enter the new input values to execute the program again, or press the BREAK key 
twice to exit from the program. 


EXAMPLE 


GAMMA TABLED 


Ås BRIGHT TAIL PROB 56,2,.1 
A=S6 B=2 RT TAIL PROB=8.1 GAMMA TABLED=69. 918511 


A,B, RIGHT TAIL PROB 12,2.35,.8 
A212 B=2.35 RT TAIL PROB=0.8 GAMMA TABLED=7. 4618437 


As Bs RIGHT TAIL PROB 6.23,3.21,.75 
Az6.23 B=3.21 RT TAIL PROB=6.75 GAMMA TABLED=2.9466713 


AB RIGHT TAIL PROB 
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TABLED VALUES 
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NORMAL TABLED 


DESCRIPTION 


This program will calculate an approximate tabled value for the normal distribution. The 
right tail probability is required input. Tabled value and the probability are displayed as 
output. 


SE 23 «= 53 


The tabled value is calculated with a maximum absolute error of 1.4 X 10 © for the right 
tail probabilities between .01 and .99. For right tail probabilities between .0001 and .01, 
or between .99 and .9999, maximum error is 5 X 10 5. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


METHODS 


An approximate formula (Reference 1, equation 26.2.23) is used to obtain an initial esti- 
mate. A one step iteration method based upon Reference 2 is then used to improve the 
accuracy. 


Dub oo IE Sd cs sa 


) 


REFERENCES 


1. Abramowitz, M., and Stegun, |. A., editors, Handbook of Mathematical Functions with 
Formulas, Graphs, and Mathematical Tables, National Bureau of Standards Applied 
Mathematics Series No. 55, Washington, D. C: Government Printing Office. 


2. Hill, G. W., and Davis, A. W., "Generalized Asymptotic Expansions of Cornish-Fisher 
Type", Annals of Mathematical Statistics, 39, 1968, pp. 1264—1273. 


VARIABLES 


23 


VARIABLE USAGE 


Significance level 


3 


Scratch 
Scratch 


=f 
2 
< x 4 


Tabled normal 


uc 
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TABLED VALUES 


NORMAL TABLED | 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 


keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 


1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 


THE PROGRAM NUMBER YOU WANT. 
3. Enter 19, press RETURN. 
4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
Enter FIND 20; press RETURN. 
. Enter OLD; press RETURN. 
. When the I/O light goes out, execute the program. 


RYN 


PROGRAM EXECUTION 


1. RIGHT TAIL PROBABILITY is displayed. Enter the probability value and press 


RETURN. 


2. The probability and tabled value are displayed. 


3. Other values may be entered to repeat execution of the program, or an exit may be 


called by pressing BREAK twice. 


EXAMPLE 
NORMAL TABLED 


RIGHT TAIL PROBABILITY .85 
P(X) 1.644688 >= 0.05000000 


RIGHT TAIL PROBABILITY .01 
P(X) 2.386623 >= 0.01000000 


RIGHT TAIL PROBABILITY .025 
P(X) 1.957355 >=  0.02500000 


RIGHT TAIL PROBABILITY . 000i 
P(X) 2.998601 >= 0.00815008 


2-14: RIGHT TAIL PROBABILITY 


C | 


- 


ES C 


i 


Plot 50: Stat Vol 1 | 


TABLED VALUES 


o 


BETA TABLED 


3 


DESCRIPTION 


This program calculates the tabled value for a beta distribution and is the inverse of the 
Beta Distribution program in the DISTRIBUTIONS section. 


a) 


Parameters a and b and the right tail probability are required input values. 


ee alium À 


The tabled value is calculated to a minimum 5 digit accuracy. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


PROGRAM LIMITS 
Parameters a and b must be greater than 0.01 for all values of P between .0005 and .9995. 


) 


METHODS 


An initial estimate for the tabled value is calculated. A quasi-linear interpolation scheme 
then uses the estimate to calculate the tabled value. 


REFERENCES 


1. Aroian, T. A., "Continued Fractions for the Incomplete Beta Function", Annals of 
Mathematical Statistics, 12 (1941), pp. 218—223. 


2. Kopitzke, R. W., “Unpublished Notes.” 


3. Peizer, D. B., Pratt, J. W., "A Normal Approximation for Binomial, F, Beta, and Other 
Common Related Tail Probabilities, |” Journal of the American Statistical Association, 
63 (1968), pp. 1457—1483. 
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BETA TABLED 


VARIABLES i 
VARIABLE USAGE 
Scratch | 
A1 Scratch 
A2 Scratch |l 
B Scratch 
B1 Scratch 
B2 Scratch i 
F1 Parameter a 
F2 Parameter b |! 
| Counter 
P Right tail probability i 
PO Scratch 
P1 Scratch a E 
Q Tabled beta V = 
T Scratch | . 
T7 Flag I| 
T8 Flag 
U Scratch M 
V Scratch 
V1 Scratch | 
X Scratch 
X1 Approximate root L 
x2 Approximate root 
X3 Approximate root | 
Y Scratch U 
Y1 Scratch . 
Y2 Scratch i 
Y3 Scratch 
U 
A 
2 7 1 6 Piot 50: Stat Vol 1 
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BETA TABLED 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


mo gg X a 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 20; press RETURN. 
4. When the I/O light goes out, execute the program. 


a a a 


3 


MANUAL LOADING FROM THE KEYBOARD 

. Insert the program tape. 

. Enter FIND 21; press RETURN. 

. Enter OLD; press RETURN. 

. When the I/O light goes out, execute the program. 


) 
RON 


een 


PROGRAM EXECUTION 


1. A, B, RT TAIL PROBABILITY is displayed. Enter the values, using commas to sepa- 
rate the entries. 


2. Input and tabled beta values are displayed. 


3. The program may be executed repeatedly by entering other values, or an exit may be 
called by pressing the BREAK key twice. 


EXAMPLE 


BETA TABLED 


AL Ba RT TAIL PROBABILITY 5,63,2.48,.6 
A=5.63 B=2.48 RT TAIL PROB=0.6 TABLED BETA=0.2490291 


A,B RT TAIL PROBABILITY 59,1.56,.3 
$259 Bai.56 RT TAIL PROB=0.3 TABLED BETA=0.031645 


A B RT TAIL PROBABILITY 18,28,.01 
ART PAÍS RT TATL PROBLO.O1 TABLED BETA=O. 7652199 


As Ba RT TAIL PROBABILITY 


Plot 50: Stat Vol 1 2 -1 7 


a 23 23 
2 


-— 
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R BY C CONTINGENCY TABLED 
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DESCRIPTION 


This program performs a contingency table analysis on a two way table. The null hypothesis 
tested is that the rows and columns are independent, against an alternative hypothesis of 
dependency. 


The number of rows and columns and all observed frequencies are required input values. 
All input data may be edited. 


Any or all of the following output options may be selected: 


1. 


o c PN 


HARDWARE REQUIREMENTS 


LL E 


=, TE ES ES 


List of data; 
Column numbers, totals, and probabilities; 
Row members, totals, and probabilities; 


Overall total; 


p. e 


Expected cell frequencies; and 

Pearsons coefficient of contingency, number of cells with expected number less than 
5, number of cells with expected number less than 1, and chi-square test with degrees 
of freedom and significance level. 


C 


A Tektronix 4050-Series Graphic System is required. If the Data Communications Interface 


(Option 1) is installed, minimum acceptable memory configuration is 16K bytes. 


METHODS 


The formulae below are applied for a contingency table with c columns, r rows, and n 


Cc E. 


observations: 


= observed frequency in row i, column j. 


[— 


nj = total of column j = nj tnj ERE 

n, = total of rowi=n,, tn, *... tni. i 
e; = expected number of row i, column j=(n, nn 

chi-square = X = (nij — eij)? /eij | 
coefficient of contingency = LI 


degrees of freedom = (r — 1) (c — 1). 


i 


REV. A, FEB. 1977 Plot 50: Stat Vol 1 ) i 


TABLED VALUES 


R BY CCONTINGENCY TABLED 


REFERENCES 


1. Mood, A., Graybill, F., Boes, D., /ntroduction to the Theory of Statistics, New York: 
McGraw-Hill, 1974. 


2. Siegel, S., Non-Parametric Statistics for the Behavioral Sciences, New York: McGraw- 
Hill, 1966. 


3. Kendall, M. G., and Stuart, A., The Advanced Theory of Statistics, New York: Hafner 
Publishing Company, 1973. 


VARIABLES 


VARIABLE USAGE 


Number of columns in the table 
Scratch 

Actual number of columns in the table 
with allowance for totals. (C+1) 


G2 Scratch 

C3 P(X? ,DF) 

C4 Scratch 

C5 Number of cells with expected value 
less than 5. 

C6 Chi-square 

C7 Number of cells with expected value 
less than 1. 

C8 Scratch 

C$ Scratch 

D The contingency table (F1,C1) 

D1 Scratch 

D2 Scratch 

D3 Scratch 

D4 Scratch 

D5 Scratch 

D6 Coefficient of contingency 

D7 Scratch 

D8 Degrees of freedom 

D9 Scratch 


Number of rows in the table 
Actual number of rows in the table with 
allowance for totals. (f+1) 
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FUNCTION KEYS 


CONTINGENCY TABLES TAPE 1 


SHIFT KEYS 
ENTER CORRECT DISPLAY ROW COLUMN 
DATA DATA DATA STATISTICS |) STATISTICS 
SHIFT KEYS 


EXPECTED 


KEY FUNCTION 


ENTER DATA allows entry of data from keyboard. 
CORRECT DATA permits correction of incorrectly entered data. 
DISPLAY DATA displays data which has been entered. 

ROW STATISTICS | calculates and displays row statistics. 
COLUMN calculates and displays column statistics. 


STATISTICS 


EXPECTED calculates and prints the expected number of observations 


VALUES for each cell. 


TESTS calculated and displays chi-square test of independence, de- 
grees of freedom, significance level, Pearsons coefficient 
of contingency, and the number of cells with expected 


values less than 5 and less than 1. 
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TABLED VALUES 


R BY CCONTINGENCY TABLED 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 


1. 
2. 


Insert the program tape. 


Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


. Enter 21; press RETURN. 


4. When the 1/0 light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 


RUN 


Insert the program tape. 
Enter FIND 22: press RETURN. 


. Enter OLD; press RETURN. 
. When the 1/0 light goes out, execute the program. 


PROGRAM EXECUTION 


1. 


Press ENTER DATA (key No. 1). ENTER THE NUMBER OF ROWS AND COLUMNS 
OF THE TABLE is displayed. Enter the number of rows, a comma, and the number 
of columns. Press RETURN. 


. ENTER THE NUMBER OF OBSERVATIONS FOR EACH CELL is displayed, and 


FOR ROW i, COLUMN j appears on the screen. Enter the observation and press 
RETURN for each cell entry. 


. The entry procedure continues until ENTRY OF TABLE COMPLETE is displayed 


to indicate that all data has been entered. 


If errors which cannot be corrected with the CLEAR key are made during data entry, 
continue to enter the data. After all data has been entered, use CORRECT DATA 
(key No. 2) to change the incorrect values. 
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CORRECT DATA (Key No. 2) Lj 


Key No. 2 allows correction of incorrectly entered values, but must not be used until the 
table has been filled. To correct an entry in the table: 


1. Press CORRECT DATA. 
CORRECT THE TABLE: ENTER ROW AND COLUMN OF CELL 
ROW, COLUMN is displayed. 


2. Enter the incorrect row number, a comma, and the incorrect column number. Press 
RETURN. 


3. The value of the existing cell entry is shown and ENTER NEW VALUE is displayed. 
4. Enter the correct value and press RETURN. 
This correction procedure continues until another function key is pressed. When all cell 


values are correct, function keys 3, 4, 5, 6, and 7 may be used to perform analysis of the 
table. 


[— EF. G G G 


C 


DISPLAY DATA, ROW and COLUMN STATISTICS, EXPECTED VALUES, and TESTS 
function key usage is demonstrated in the examples which follow. 


E: 


EXAMPLE 1 


Calculate the chi-square test for the following table: 


Column Column Column 


C O E X 


Load the program and press ENTER DATA (key No. 1). All data for row 1 is entered, then 
the data for row 2, and row 3. In this example, the value for row 2, column 1 should have 
been 7, but was entered as 15 in order to illustrate the data correction procedure. 


C 
oL. Es 
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The table is entered as shown below. 


Rx € Contingency Table Analysis 


Enter the number of rows and columns of the table 3,3 


Enter the number of observations for each cell 


for row 1, column 1: 4 
for row 1, column 2º 5 
for row 1, column 3: 6 
for row 2; column i: 15 
for row 2, column 2: 6 
for row 2, column 3: 9 
for row 3, column 1: 12 
for row 3, column 2: 13 
for row 3, column 3: 14 


Entry of table complete 


d a ek se sen ed 


2 


To correct the entry in row 2 column 1, press CORRECT DATA (key No. 2) 


Correct the table: Enter row and column of the cell 
Row, column 2,1 

The existing cell entry is 15 

Enter neu value ? 


Rows column 


dr 23 


Press DISPLAY DATA (key No. 3) and the data is displayed as follows: 


d LL 23 


— 
— 


Data Table 


20 
-0 

£ 

64 
= 
"m 


14 


2y d coc 
så 
€ 
nm 
o 
o 
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Row statistics are displayed when ROW STATISTICS (key No. 4) is pressed. E 


Row Statistics 


Row N P 
1 15 6. 197368421853 i 
2 22 0,289473684211 
3 39 9.513157894737 I 
Total 76 
I 
Column statistics are displayed when COLUMN STATISTICS (key No. 5) is pressed. LU 
Column Statistics | 
| 
NTG 
Colunn N P 
1 23 0. 302631578947 å 
2 24 0.315789473684 
3 29 A. 381578947368 U 
Total 76 
Press EXPECTED VALUES (key No. 6) and the expected number of observations in each 
cell is printed. H 
Expec ted Values i 
Row 1 
4.53947368421 4.?3684210526 5. 72368421053 i 


Row 2 i 
6.6535789473684 6.947368421085 8.739473684211 \ j 


Rou 3 
11.8026315789 12.3157894737 14.8815789474 
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Press TESTS (key No. 7) and the Chi-Square value, degrees of freedom, significance level, 
and the number of cells with small observed or expected frequencies are printed. 


Tests 


Chi Square 
DF 
P<Chi-Square,s DF) 


9.376811994003 
0.98439 


Coefficient of Contingency 8,607016529576935 
# of cells with exp. val, « 8 2. 
* of cells with exp. val. € 1 G. 


nu 


The chi-square test of .376 with significance level .98 indicates that the hypothesis of row 
and column independence should not be rejected. The table shows that the number of ob- 
servations in each row increases in relation to column number increase, and that an ob- 
servation in a given row provides no specific information about column placement of the 
observation. 


EXAMPLE 2 


School children were classified according to standard of clothing and according to intelligence. 
Standards were as follows: 


ED oJ en | 


= | -3 
) 


CLOTHING STANDARDS INTELLIGENCE STANDARDS 


f A very well clad A slow and dull 
B well clad B dull 
n C poor but passable C slow but intelligent 
` D very poorly clad D fairly intelligent 
E distinctly capable 
n F very able 
1 Results of the classification are shown below: 


INTELLIGENCE 
CLOTHING | A D 


D 


— LI. 22 
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The program was loaded and the data entered in the following manner: 


Enter the number of rows and columns of the table 4,6 


Enter the number of observations for each cell 


for rou 
for row 
for row 
for rou 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 
for row 


column 
column 
column 
column 
column 
column 
column 
column 
column 
column 
column 
column 
column 
column 
column 
coluan 
column 
column 
column 
column 
column 
column 
column 
column 


Rx C Contingency Table Analysis 


Av Oe TO GU 49 9» ^P o9 na O9 P9 Un Da qu 99 nu na an unas "D vo os o 


ON CDI da O4 P0 > OM CID GID OM OG GJ P3 OO & GJ PO ra 


Entry of table complete 


To display the table, press DISPLAY DATA (key No. 3). 


Data Table 


48 
39 


?0 


22 
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Press ROW STATISTICS (key No. 4) 


Row Statistics 


TABLED VALUES 


R BY C CONTINGENCY TABLED 


Row N P 
1 636 8. 368695652174 
2 751 9.435362318841 
3 265 6. 153623188466 
4 73 6. 0423132455797 
Total 1725 


COLUMN STATISTICS (key No. 5) produces the result shown below: 


Column Statistics 


Column N P 
1 138 0.0753623188406 
2 219 Ö. 126956521739 
3 487 8. 235942628986 
4 535 0.310144927536 
5 375 0.2173913647345 
6 59 9.0342028985567 
Total 1725 


EXPECTED VALUES (key No. 6) will cause the display of the expected number of ob- 
servations for each cell. 


Expected Values 
se i 


dr. 9304347928 
138. 260969565 


80. 7443473261 
21. 7530434793 


150. 053120425 


Row 2 

56. 5971014493 25, 3443478261 17?. 192463768 232.91884858 
163.260869565 25.6963768116 

ov 3 Vane : 
19,97 10144928 33.6434782609 62. 5246375812 $2.1884957971 
7. 6886956522 3.06376811594 

Rou 4 " : 
5.50144927536 9,26782602696 17.2237691 159 22.6405797 181 

2.4968115942 


15.8695652174 
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TESTS (key No. 7) may be used to display the test results which follow: 


Tests 


Ec. E 


se Square = DE 
P(Chi-Square, DF > = 0.88000 
Coefficient of Contingency = 0.303347924145 
8 of cells with exp. val. € 5 = f. 

# of cells with exp, val. < 1 = 6. 


The chi-square test of independenceis 174.8 with a significance level less than 1076, so we 
reject the hypothesis that clothing standards and intelligence are independent. The row 
location of an individual provides some information about the individuals column location. 
Clothing standard A, for example, indicates that the individual is most likely to be in in- 
telligence class D or E. If his clothing standard is D, he is most likely to have intelligence 
classification A or C. 


C A GG O QG A E E. Es 


| 
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DESCRIPTION 


This program will analyze a contingency table with 2 rows and 2 columns. The four data 
points required as input values are used to generate optional output as follows: 


1. Data list, expected number of observations, cell probabilities, and table of expected 
values from observed values; and 


2. Chi-square test with and without Yates’s correction for continuity, significance level, 


and Fisher’s exact test. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


METHODS 
Let Nii i = 1,2; j = 1,2 bethe table entries. 


Let Ni. = Nii*Ni2 and 


Nj = Nij*N5j be the respective row and column totals, with N the 


total number of observations. 
The expected number in row i, column j is 
Ej = NiN j/N 
Probabilities are Ei/N- 
Deviations are Ni-Eij. 
Let D = N4.N5 N 4N 9 , and 
VENN MING: 


The chi-square values calculated are: 


chi-square uncorrected = 
N V?, and 
D 


chi-square corrected = N(|V|-N/2)? /D. 
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Fisher's test is a one tail test for positive association in a table. It is the cumulative prob- 
ability of the observed table and all less likely tables under the assumption of fixed marginal 
totals. 


T 


The chi-square tests are two-tail tests of independence (pp. 571—573, reference 3). Cochran 
(Reference 1) recommends that the continuity corrected chi-square test be used if total 
number of observations is at least 40, or if there are at least 20 observations and no cell has 
an expected frequency less than 5. 


C C 


--—" 
— 


Corrected and uncorrected chi-square tests and the Fisher's test are printed so an appropriate 
choice may be made for any data. 


C 


CZ 


REFERENCES 


1. Cochran, W. G., "Some Methods for Strengthening the Common Chi-Square Tests”, 
Biometrics, 1954, Vol. 10, pp. 417. 


E 


2. Fisher, R. A., Statistical Methods For Research Workers, 1946, Edinburgh: Oliver and 
Boyd. 


q 


3. Kendall, M. G., Stuart, A., The Advanced Theory of Statistics, Vol 2, New York: 
Hafner, 1973. 


VARIABLES 


CO O 


VARIABLE USAGE 


Scratch variable for chi-square 


E CGC 


A1 Scratch variable for chi-square 

A2 Scratch variable for chi-square 

A3 Scratch variable for Fishers test 

B Scratch variable for chi-square 

B1 Scratch variable for chi-square H 

B2 Scratch variable for chi-square 

C1 Scratch variable for chi-square | 

C9 Scratch variable for chi-square 

D Data matrix 

F Scratch i 

F1 Scratch 

| Scratch i 

K Scratch U 
Piot 50: Stat Vol 1 |) 
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M Scratch 

N Scratch 

P1 General row index 

P2 General Column index 
P6 Chi-square value 

x Scratch 

X1 Scratch 

X2 Scratch 

X3 Scratch , 
Y Chi-square probability 


FUNCTION KEYS 


CONTINGENCY TABLES TAPE 1 


SHIFT KEYS 


ENTER sea DISPLAY 
DATA DATA DATA STATISTICS STATISTICS 


SHIFT KEYS 


EXPECTED 


KEY KEY FUNCTION 


1 ENTER DATA allows data to be entered from the keyboard. 

2 CORRECT DATA | permits correction of data 

3 DISPLAY DATA 

4 ROW STATISTICS] keys 3, 4, 5, and 6 perform identical functions. Pressing 

5 COLUMN any of these keys will cause the data, expected frequencies, 
STATISTICS row statistics, and column statistics to be displayed. 

6 EXPECTED 
VALUES 

7 TESTS displays Chi-square and Fisher's tests. 
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PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 22; press RETURN. 
4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 23; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the I/O light goes out, execute the program. 


PROGRAM EXECUTION 
1. Press the ENTER DATA key. The message FOR ROW 1, COLUMN 1 is displayed. 


2. Enter the data value and press RETURN. The display screen will then request the 
next value. 


3. Repeat step 2 until all four entries have been completed and the message ENTRY OF 
TABLE COMPLETE is displayed. 


4. Press key 3, 4, 5, or 6 to display the descriptive statistics. Key No. 7 will cause the 
chi-square tests and the Fisher’s test to be displayed. Data corrections may be ac- 
complished by pressing key No. 2. 


5. To re-execute this program with other values, type RUN and repeat the program 
execution procedure. 
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EXAMPLE 1 
Perform the chi-square test for independence in the following table: 


Column 1 Column 2 


Load the program and press ENTER DATA (key No. 1). 


& X 2 CONTINGENCY TABLE 


ENTER THE NUMBER OF OBSERVATIONS FOR EACH CELL 


FOR ROM 1,COLUMN 1: 2 
FOR ROW 1,COLUMN 2: 3 
FOR ROW 2,COLUMN 1: 4 
FOR ROW 2,COLUMN 2: 2 


ENTRY OF TABLE COMPLETE 


oo 13 LL —J3 —3 3 3 —3 
D 
o 
= 
N 
q 
N 


) 


In the example, the entry in row 2, column 2 is incorrect. Å 4 was entered (instead of 5) to 
demonstrate the data correction option. 


23 


E ao 


Data correction is accomplished as shown below: 


Press CORRECT DATA (key No. 2). 


z3 


TO CORRECT AN ENTRY OF THE TABLE, ENTER ROW 
AND COLUNN OF THE CELL. 


ROW, COLUMN = 2,1 
THE EXISTING CELL ENTRY IS 4 


ENTER THE NEW VALUE: 5 
ROW, COLUMN = 


ces 33 
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The display shown below appears when DISPLAY DATA, ROW STATISTICS, COLUMN 
STATISTICS, or EXPECTED VALUES (keys 3, 4, 5, or 6) are pressed. 


DATA TABLE 
ROW 1: 2 3 9 
ROW 2: 5 2 ? 
TOTAL : ? s 12 


EXPECTED VALUE TABLE 


2. 91666666667 2 
4.08333333333 2 


666556 
PROBABILITIES 
0. 243055555556 $ 173611111111 0.41666666b667 
0.348027777?7?8 . 243055555556 8, 583333333333 
0. 583333333333 à. 41666666665? 1.8 


TABLE OF DEVIATIONS 
-0.916666666667 8 


+31 

0.91666666666? -0.9 

Test are displayed as shown here when TESTS (key No. 7) is pressed. 
CHI-SQUARE AND FISHER’S TEST 


CHI-SQUARE UNCORRECTED= 1, 18530612245 
DEGREES OF FREEDOM=1 
SIGNIFICANCE LEVEL=0, 27627744844 


CHI-SQUARE CORRECTED= 0.244897953184 
DEGREES OF FREEDGM=1 


SIGNIFICANCE LEVEL=0. 620690717672 


FISHER'S TEST=8. 318666666657 
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Both chi-square tests have high significance levels and the hypothesis of independence would 
not be rejected. 


The significance level for the corrected chi-square test is approximately twice that of the 
Fisher's exact test. This is because the chi-square corrected is a two tail test and the Fisher's 
is a one tail. 


EXAMPLE 2 
The following table classifies 42 children by feeding method and dental 
occlusion: 
Normal Teeth Mal-occluded 
Breast-fed 4 16 
Bottle-fed 1 21 


Data entry appears as shown. 


2 X 2 CONTINGENCY TABLE 


ENTER THE NUMBER OF OBSERVATIONS FOR EACH CELL 
FOR ROM i COLUMN 1: 4 
FOR ROM 1,COLUMN 2: 16 
FOR ROW 2,COLUMN 1: 1 
FOR ROK 2,COLUMN 2: 21 


ENTRY OF TABLE COMPLETE 
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TEST AND REGRESSIONS 


ONE SAMPLE t-TEST 


DESCRIPTION 


This program performs the one sample t-test of a population mean, determines confidence 
intervals, and calculates data sample statistics. The mean tested is the hypothesized value 
for a single sample of data. 


Input is single data points which may be entered from the keyboard or from magnetic tape. 
Keyboard entered data may be stored on the program tape or on data tapes, and data saved 
on tape may be deleted, corrected, or added to. 


Output may be selected from the following options: 


. list of data entered; 


. t-test of the mean equal to a hypothesized value; 


1 
2 
3. confidence intervals on the mean for 90°, 95°, and 99° levels of confidence; 
4. histogram of data; and 

5 


. descriptive statistics. 


Descriptive statistics which may be displayed are: 


. number of observations; 
. sample mean, variance, and standard deviation; 


. standard error of the mean; 


. coefficient of variation; 


1 
2 
3 
4, data minimum, maximum, and range; 
5 
6. skewness; and 

7 


. kurtosis. 


HARDWARE REQUIREMENTS 
A Tektronix 4050-Series Graphic System is required. 


The number of data points which may be stored internally is dependent upon the memory 


size of the machine being used. The following table shows the maximum observations per- 
mitted in internal storage: 
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MEMORY SIZE MAXIMUM NUMBER OF 
OBSERVATIONS 
8K 70 
16K 400 
24K 1400 
32K 2400 


The number of data points in each file should not exceed the appropriate limitation defined 
in the previous table. 


If machine capability is enhanced through the addition of additional optional memory, vari- 
able WO should be reset. This variable establishes the data matrix dimension and is located 
at line 110 of the file 28 program listing. If the value of WO is altered, data files must be 
updated to reflect the increased capability. 


PROGRAM LIMITS 
The number of cells specified for histograms should not exceed 20. 


C 


TAPE STRUCTURE 


All data tapes used by this program must be binary coded, and each must contain a Header 
File and one or more data files. The Header File, used to load data into memory, includes 
the file number assigned to the first data file, and the number of files in the data set. 


Each data file contains a figure showing the number of data points in the file. This number 
is followed by data points and calculated data for each point. 


The following illustrates the one sample t-test file structure on tape. 


HEADER FILE DATA FILE | DATA FILE 


E Number of ` Number Data | Data Data Numbers 
ii pam Files in of Point | Point Point of 
Data Set Points Xi x, X, Points 


1 
3 -2 ! Plot 50: Stat Vol 1 


nme esse ene Vane | 
) 


3 = 


LB ^ OG 35 02 23 
2 


2 


SS 


am oJ 3 223 


TEST AND REGRESSIONS 
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METHODS 


Formulae used in this program are conventional and are shown in the Appendix. 


REFERENCES 
1. Dixon, W. J., Massey, F. J., Introduction to Statistical Analysis, New York: McGraw- 
Hill, 1957. 


2. Dunn, O. J., Clark, V. A., Applied Statistics: Analysis of Variance and Regression, 
New York: John Wiley and Sons, 1974. 


3. Peizer, D. B., Pratt, S. W., "A Normal Approximation for the Binomial, F, Beta, and 
Other Common Related Tail Probabilities |, Journal of The American Statistical 


Association, 1968. 


4. Snedecor, G. W., Cochran, W. G., Statistical Methods, sixth edition, Ames, lowa: lowa 
State University Press, 1967. 
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VARIABLES 


VARIABLE USAGE 


Data matrix 


Calculate significance levels 
Calculate significance levels 
Calculate significance levels 
Calculate significance levels 
Calculate significance levels 

Save summary statistics 
Maximum for histogram generation 
Cell width for histogram 

Scratch for histogram 

Scratch for histogram 

Scratch for histogram 

Scratch for histogram 

Scratch for histogram 

Scratch for histogram 

Scratch for histogram 

Matrix for storing histogram data 
Flag for saving data on tape 


O = don't save 
1 = data on program tape 
—1 = data on data tape 
Flag indicating tape 
0 = data tape in 
1 program tape in 
File required for overlaying program 


Stop entry performed 

Program file currently in memory 
Overlay key pressed 

String variable for entering data 
Beginning file number and number 
of files used 

Scratch 

Scratch 

Scratch 

Scratch 


3-4 
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VARIABLE USAGE 


Number of observations in current 


file 
N1 Scratch 
N2 Scratch 
P Scratch 
Q Scratch 
S Scratch 
T Scratch 
WO Dimension of A, the data matrix 
X Scratch 
X1 Scratch 
Scratch 
X3 Scratch 
X4 Scratch 
^ Y1 Scratch 
Y4 Scratch 


FUNCTION KEYS 


SHIFT KEYS 


DATA STOP 
FROM TAPE ENTRY 


ENTER ADD DELETE CHANGE LIST å 
DATA DATA POINT å POINT 


ILL DM M M NM M He HE HE Ee ==" M 
x 
N 


É EH Em | 
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KEY FUNCTION 

1 ENTER DATA allows entry of data from keyboard. 

2 ADD DATA allows data to be added to an existing file. 

3 DELETE POINT permits data point to be deleted from memory and tape. 

4 CHANGE POINT allows data value to be changed. 

5 LIST DATA lists data previously entered. 

6 HISTO causes histogram of data to be displayed. 

7 STATISTICS calculates and displays the number of observations, the 

sample standard deviation, standard error of the mean, data 

minimum, data maximum, data range, coefficient of var: 
ation, skewness, and kurtosis. 

8 T-TEST CI calculates t-test and confidence intervals for the hypo- 
thetical mean. 

11 DATA FROM TAPE | allows entry of data from taped file. U 

12 STOP ENTRY terminates data entry from keyboard and permits storage 


of previously entered data. 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 
1. Insert the program tape. 
T 


— m 
-H 


2. Press AUTO LOAD. The program directory will be displayed with the message EN 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 23; press RETURN. 


4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 24; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the I/O light goes out, execute the program. 
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PROGRAM EXECUTION 

Program execution procedures are determined by the input method and data storage needs. 
The 5 options from which a selection must be made are as follows: 

Input data from the program tape; 

Input data from a previously created data tape; 

Input from the keyboard, data to be stored on the program tape; 


Input from the keyboard, data to be stored on a data tape; and 


a co 19: 3 


Input from the keyboard, no data to be stored on tape. 
Execution procedures are described for each option. 


OPTION 1, EXECUTION WITH DATA INPUT FROM THE PROGRAM TAPE 

Execute the program using a data set previously stored on the program tape. 
1. Press DATA FROM TAPE (key No. 11). IS DATA ON PROGRAM TAPE is displayed. 
2. Enter YES from the keyboard; press RETURN. When all data from the program tape 


has been read into memory, the I/O light goes out and ALL DATA ENTERED is 
displayed. 


3. Use the appropriate function keys to execute program operations. 


OPTION 2, EXECUTION WITH DATA INPUT FROM DATA TAPE 


To execute the program with input from an existing data tape: 


1. Press DATA FROM TAPE (key No. 11). IS DATA ON PROGRAM TAPE is displayed. 
2. Enter NO from the keyboard; press RETURN. HEADER FILE is displayed. 


3. Enter the numeric designation of the header file. Press RETURN. The alarm bell sounds 
and INSERT DATA TAPE, PRESS RETURN is displayed. 


4. Insert the correct data tape and press RETURN. When the data transfer is complete, 
the 1/0 light goes out and ALL DATA ENTERED is displayed. 


5. Use the appropriate function keys to execute program operations. 


OPTION 3, EXECUTION WITH KEYBOARD INPUT, DATA TO BE STORED ON 
PROGRAM TAPE 


To execute the program and store keyboard entered data on the program tape: 


1. Press the ENTER DATA (key No. 1). DO YOU WISH TO STORE DATA ON TAPE 
is displayed. 
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. Enter YES from the keyboard; press RETURN. DO YOU WISH TO STORE DATA 


ON PROGRAM TAPE is displayed. 


. Enter YES from the keyboard and press RETURN. 


ENTER DATA 
X = is displayed. 


. Enter the value of X and press RETURN for each data point. 


5. When all data has been entered, use the appropriate function keys to perform other 


OPTION 4, EXECUTION WITH KEYBOARD INPUT, DATA TO BE STORED ON A DATA 


program operations. Data storage will be automatic. 


. If no analysis is desired, press the STOP ENTRY (key No. 12) to store the data and 


update the Header File. 


TAPE 


To execute the program and store keyboard entered data on a data tape, follow the steps 
below: 


T. 


Press the ENTER DATA (key No. 1). DO YOU WISH TO STORE DATA ON TAPE 

is displayed. 

Enter YES from the keyboard and press RETURN. DO YOU WISH TO STORE DATA 
ON PROGRAM TAPE is displayed. 


. Enter NO from the keyboard and press RETURN. HEADER FILE is displayed. 


4. Enter the number of the header file and press RETURN. INSERT DATA TAPE, 


PRESS RETURN is displayed and the alarm bell sounds. 


. Insert the correct data tape and press RETURN. 


ENTER DATA 
X = is displayed. 


6. Enter the X value and press RETURN for each data point. 


7. When all data points have been entered, use the appropriate function keys to perform 


program operations. Data storage wil! be automatic. 


. If no data analysis is desired, press the STOP ENTRY (key No. 12) to store the data 


and update the Header File. 
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OPTION 5, EXECUTION WITH KEYBOARD INPUT AND NO DATA STORAGE 
To execute the program without storing data for later use, proceed as follows: 
1. Press ENTER DATA (key No. 1). DO YOU WISH TO STORE DATA ON TAPE is 
displayed. 
2. Enter NO from the keyboard and press RETURN. 
ENTER DATA 
X = is displayed. 
3. Enter the X value and press RETURN for each data point. 


4. When all data points have been entered, use the function keys to perform any of the 
operations described on the next pages. 


PROGRAM OPERATIONS 


Upon completion of the initial data entries, operations are performed with the one sample 
t-test function keys. The operations described here may be performed in any sequence. 


ADD DATA 


ADD DATA (key No. 2) will permit addition of new data to the file currently in memory. If 
the file also resides on tape, the alarm bell sounds to indicate that the tape should be inserted. 
A message on the screen will identify the kind of tape required; program tape or data tape. 


When the tape has been inserted, press RETURN. 


ENTER DATA 
X = will be displayed, and values may be entered. The RETURN key must be 
pressed after each value has been keyed. 


Other operations may now be performed, or the STOP ENTRY (key No. 12) may be used 
to update the file. 


DELETE POINT 


DELETE POINT (key No. 3) allows a value to be deleted from an existing file. The data 
point will be deleted from internal storage and from the tape file. All subsequent calculations 
will exclude the deleted point. 


If the alarm bell sounds when DELETE POINT is pressed, insert the tape indicated by the 
message on the screen, and press RETURN. 


When DELETE VALUE appears on the screen, enter the value to be deleted and press RE- 
TURN. The deletion will be completed with no further message. If the value entered is not 
in the stored data, DATA NOT FOUND is displayed and no deletion occurs. 3-9 
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CHANGE POINT 


Incorrect data points are corrected with CHANGE POINT (key No. 4). When the key is 
pressed, the alarm bell and a message on the screen indicate any tape requirement. If the 
proper tape has been inserted, the screen will display CHANGE: INCORRECT VALUE, 
CORRECT VALUE. When the incorrect value, a comma, and the correct value are entered 
and the RETURN key is pressed, correction will occur immediately. 


If the value specified is not in the data file, DATA NOT FOUND will appear on the screen. 


LIST DATA 


To list data previously stored, press LIST DATA (key No. 5). If a tape insertion is needed, 
the alarm bell will sound and a message will be displayed. Insert the correct tape and press 
RETURN. 


DATA IS NOT AVAILABLE TO BE LISTED will be displayed if the data is not stored. 
Otherwise, the data will be listed on the screen. 


C 


|f the screen becomes full, the listing will be continued by pressing the PAGE key. 


HISTO 


Stored data points are used to generate the histogram displayed by pressing HISTO (key No. 
6). The alarm bell and a displayed message indicate any necessary magnetic tape insertion. 
RETURN must be pressed after the correct tape has been inserted. 


When HISTO (key No. 6) is pressed, the data minimum, maximum, and range are displayed 
with the following message: 


SELECT OPTION 


1 SPECIFY FOF CELLS 

2 SPECIFY CELL WIDTH 

3 SPECIFY OFFSET AND CELL WI DTH 
ENTER OPTI ON 


After entering the option number, press RETURN. The screen will then display either # OF 
CELLS; CELL WIDTH; or OFFSET, CELL WIDTH. Enter the correct information and 
press RETURN. 


C 


The number of cells for option 1 must not exceed 20. 
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DO YOU WANT NORMAL CURVE will now be displayed. If YES is entered, a normal 
curve overlay will appear in addition to the histogram. If NO is entered, the histogram will 
be drawn without the curve. 


STATISTICS 


Descriptive statistics are displayed when STATISTICS (key No. 7) is pressed. Any required 
tape insertion is indicated by the alarm bell and a message on the display screen. 


Statistics displayed are: 
Number of observations; 
Sample mean, variance, and standard deviation; 
Standard error of the mean; 
Data minimum, maximum, and range; 
Coefficient of variation; 
Skewness; and 
Kurtosis. 


t-TEST AND CONFIDENCE INTERVAL 


t-TEST CI (key No. 8) initiates the calculation and printing of the hypothesized mean t-test 
and confidence intervals. If the alarm bell sounds when the key is pressed, insert the tape 
indicated by the display screen. When the RETURN key is then pressed, the display will be 
HYPOTHESIZED MEAN. 


After the hypothesized mean is entered and the RETURN key is pressed, the t-test and con- 
fidence intervals are displayed. 


STOP ENTRY 


STOP ENTRY (key No. 12) is only used if data has been entered and saved on tape and if 
no other analysis has been done. This key will ensure that all data is saved on tape. 


EXAMPLE 1 


Test the hypothesis that the data below is a sample from a population with mean greater 
than or equal to 29. 


27.2, 34.1, 33.9, 28.6, 29.2, 30.1, 31.2, 30.5, 33.6, 30.9, 31.3, 31.0, 32.3, 31.5, 
26.2, 29.9, 34.0, 32.1, and 30.7. 
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Load the program and begin program execution by pressing ENTER DATA (key No. 1). With 
so few observations, all data may be stored internally and there is no need to store data on 


tape. 


ONE SAMPLE T-TEST 
DO YOU WISH TO STORE DATA ON TAPE? NO 


Data entry is shown below. 


NTER DATA 
27.2 


som 


e 
^ ^ DON AG G4 OO OY OA ra 


Gå 04 GJ PI PD GI GIGI CU GI GJ GI GI 
ON BOK 7 — OG OG 0D 9 CD 


DE DC XE 3€ 3€ DC 3€ 3€ DE 5€ DE 3€ DIE X 2€ 25€ II 


Press LIST DATA (key No. 5) 


27. 2080 
38.1982 
31.3080 
34.0080 


DATA 
34.1000 
31.2000 
31.0800 
32.1008 


33. 9000 
30. 5000 
31.5008 
38. 7800 


C 


28.680008 29.2080 
33, 6808 39. 9000 
26.2000 29, 3000 


c 
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Press HISTOGRAM (key No. 6) 


HISTOGRAM 


DATA MIN = 2 
DATA MAX = 3 
DATA RANGE = 


SELECT OPTION 
1 SPECIFY # OF CELLS 
2 SPECIFY CELL WIDTH 
3 SPECIFY OFFSET AND CELL WIDTH 


ENTER OPTION! 
# OF CELLS 16 


DO YOU WANT NORMAL CURVE YES 


To specify the number of cells, enter 1, and press RETURN. In the example above, 10 


3 


S \ cells and a normal overlay were specified. 


N x REL FREQ 


> 
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Because the cell limits are confusing, a better histogram might be generated by setting the j 
offset at 26.0 and the cell width at 1. This will produce fewer than 20 cells and clear cell U 
limits. 
HISTOGRAM | 
DATA MIN = 26.2 
DATA MAX = 34.1 
DATA RANGE = 7.9 i 
SELECT OPTION ; 
1 SPECIFY 4 OF CELLS Li 
2 SPECIFY CELL WIDTH 
3 SPECIFY OFFSET AND CELL WIDTH 
ENTER OPTION3 i 
OFFSET, CELL WIDTH 26,1 
DO YOU WANT NORMAL CURVE YES I 
N x REL FREQ | 
6. 33.3 Ny L| 
; U 
4 
I 
L 
l 
f i 
+2. 600 +2.800 +3.000 +3.200 +3.400 ør 
41091 | 
3-1 4 Plot 50: Stat Vol 1 |] 
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n ONE SAMPLE t-TEST 
n When STATISTICS (key No. 7) is pressed, the following display appears: 
DESCRIPTIVE STATISTICS 
a i: i 
MEAN = 30.8888888889 
f VARIANCE = 4.94104575165 
| STD DEV = 2. 2228463176 
n DATA MIN = 26.2 
DATA MAX = 34.1 
n DATA RANGE = 7.9 
STANDARD 
ERR OF MEAN = 0.5239239091569 
n COEFFICIENT 
l OF VARIATION = 7. 19626505697 
-ġ. 494061927 
n AN SKEWNESS = 9.325 
ni KURTOSIS s 2.66356133199 
n T-TEST (key No. 8) initiates a request for the hypothesised mean. 
1 T-TEST 
HYPOTHESISED MEAN =29 
T TEST OF MEAN = 29 IS 3,60523208015 
1 WITH 17 DEGREES OF FREEDOM 
| SIGNIFICANCE LEVELS : 
ONE TAIL TEST 20.600148 
f TWO TAIL TEST =0.46296 
SAMPLE MEAN =30. 8988888589 
N CONFIDENCE INTERVALS FOR POPULATION MEAN : 
LEVEL LOWER LIMIT UPPER LIMIT 
904 i29.9772421149 , 31.6605356629) 
95% (29.7832214569 , 31.9945563209) 
99% (29,3741316463 , 32.4076461315) 
: 1 
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Because the sample mean is 30.9 and the significance level for a one tail test is .0007, we 
reject the hypothesis that the mean is less than or equal to 29. The hypothesis that the 
mean is greater than or equal to 29 is accepted. 


If the test had been two tailed, we would have rejected the hypothesis that the mean is 
equal to 29 at the .0014 significance level. 
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DESCRIPTION 


This program tests the hypothesis that the means of two populations are equal. Data 
samples from the populations may be entered from the keyboard or from magnetic tape, 
and the input may be corrected, deleted, or updated. 


Output may be selected from the following options: 


1. data list; 


2. number of observations, mean, standard deviation and error of the mean, data mini- 
mum, maximum and range, skewness, and kurtosis; 


. pooled variances; 
. t-test of equality of the means; 


. histogram for either or both samples; 


a ao fF Q 


. confidence intervals for each sample mean, and for the difference of the means. 


HARDWARE REQUIREMENTS 


A Tektronix 4050-Series Graphic System is required. If the Data Communications Interface 
(Option 1) is installed, minimum acceptable memory configuration is 24K bytes. 


The maximum number of data observations which may be stored is dependent upon the 
memory size of the machine being used. Capabilities are: 


Memory Maximum Number 
Size of Observations 


WO, the variable located in line 110 of program file 34, is set to the correct value of the 
maximum number of observations. If memory size is increased by the addition of optional 
memory modules, this variable must be corrected to conform to the table shown above. Be- 
cause the allowable number of observations per file is determined by the storage capability 
of the machine, this change will also necessitate data file updating. 
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TAPE STRUCTURE 


This program recognizes two types of taped files, both of which must be in binary format. 


The Header File consists of two numbers; the Header File designation number, and the 
number of files used for data storage. Header File information is required for loading data 
into memory. 


Data Files each contain four types of information in the order specified below: 


1. The number of observations (N,) from the first population sample; 
2. The number of observations (N3) from the second population sample; 
3. Observation data from the first population; and 


4. Observation data from the second population. 


The size of each Data File may be stated as N4 +N5+2. 


METHODS 


All input data entered from the keyboard is saved in the one dimensional array A. When the 
number of observations in A is equal to the value of WO, the entire array is either written to 
tape or erased from memory, depending upon the alternative selected during program 
execution. 


The standard formulae used in this program are contained in the Appendix. 


REFERENCES 
1. Dixon, W. J., Massey, F. J., Introduction To Statistical Analysis, New York: McGraw- 
Hill, 1957. 


2. Dunn, O. J., Clark, V. A., Applied Statistics: Analysis of Variance and Regression, 
New York: John Wiley and Sons, 1974. 


3. Peizer, D. B., Pratt, J. W., “A Normal Approximation for the Binomial, F, Beta, and 
Other Common Related Tail Probabilities |”, Journal of the American Statistical 


Association, 63, 1968. 


4. Snedecor, G. W., Cochran, W. G., Statistical Methods, Sixth Edition, Ames lowa: lowa 
State University Press, 1967. 
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USAGE 


Data matrix 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Save summary statistics 

Cell width for histograms 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

String variable for alpha information 
Matrix for storing histogram data 
Flag for saving data on tape 

O = do not save 

1 7 save on data tape 

—1 = save on separate tape 

Flag for tape type 

O = data tape in 

1 = program tape in 

File required for overlaying program 
Flag for step entry 

O = not performed 

1 = step entry performed 

program file currently in memory 
overlay key pressed 

Beginning file number, number of files used 
Scratch 

Scratch 

Scratch 

Scratch 

String variable for input alpha information 
Number of observations in the current file 
Scratch 

Scratch 
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VARIABLE USAGE 


Scratch 
Scratch 
Scratch 
Scratch 
Scratch 
Scratch 
Dimension of the data matrix 
Scratch 
Scratch 
Scratch 
Scratch 
Scratch 
Scratch 
Scratch 
Scratch 
Y4 Scratch 


FUNCTION KEYS 


T-TESTS TAPE 1 


SHIFT KEYS 


DATA 
FROM TAPE 


ENTER ADD DELETE CHANGE LIST 
DATA DATA POINT POINT DATA 
- SHIFT KEYS 


STOP 
ENTRY 
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KEY KEY 


1 ENTER DATA 

2 ADD DATA 

3 DELETE POINT 
4 CHANGE POINT 
5 LIST DATA 
6 HISTO 

7 STATISTICS 


T-TEST CI 


-3 pum | 
2 
œ 


12 STOP ENTRY 
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FUNCTION 


allows entry of data from keyboard. 

permits data to be added to an existing file. 

allows data point to be deleted from memory and tape. 
allows data value to be changed. 

lists data previously entered. 

causes histogram of data to be displayed. 


calculates and displays the number of observations, the 
sample mean, sample variance, sample standard deviation, 
standard error of the mean, data minimum, data maxi- 
mum, data range, coefficient of variation, skewness, and 
kurtosis. 


calculates t-test and confidence intervals for the hypo- 
thetical mean. 


allows entry of data from tape. 


terminates data entry from keyboard and permits storage 


| of data, if tape storage option was elected. 


This program may be loaded automatically from the directory, or manually from the 


AUTOMATIC LOADING FROM THE DIRECTORY 


z3 


1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 


pe a 


2 
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THE PROGRAM NUMBER YOU WANT. 
3. Enter 24, press RETURN. 
4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 25; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the 1/0 light goes out, execute the program. 3-2 1 
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PROGRAM EXECUTION 
If input data is to be entered from an existing magnetic tape file, follow the instructions 


concerning the DATA FROM TAPE (key No. 11). Otherwise, data must be entered from the 


keyboard with the ENTER DATA (key No. 1). 


ENTER DATA 


1. 


Press ENTER DATA (key No. 1) to display DO YOU WISH TO STORE DATA ON 
TAPE? Reply YES or NO through the keyboard and press RETURN. If you replied 
NO, omit steps 2, 3, and 4; go to step number 5 below. If you replied YES, continue 
the sequence. 


. DO YOU WISH TO STORE DATA ON PROGRAM TAPE? is displayed. If the file is 


on the program tape, type YES; press RETURN, and go to step 4. If the data is on a 
separate data tape, reply NO and press RETURN. 


. To establish the identity of the input data file, the message HEADER FILE # 


is displayed. Enter the correct Header File number and press RETURN. 


. If the alarm bell sounds, insert the tape indicated by the message on the screen and 


press RETURN. 


. The display screen shows: 


ENTER DATA 
X SAMPLE 


Enter the X value and press RETURN. Repeat this step for each value in the first 
sample. To proceed to data entry for the second sample, press RETURN without 
entering a data value. 


. The message 


Y SAMPLE 


will be displayed. Enter the Y values and press RETURN after each entry. When all 
Y values have been entered, use the function keys to execute data analysis. 


DATA FROM TAPE 


1. To allow entry of data from a previously created magnetic tape file, press DATA FROM 


TAPE (key No. 11). 


Data tapes used by this program must conform to the tape standards explained 
in the TAPE STRUCTURE section. 
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TWO SAMPLE t-TEST 


If the program tape has been removed from the machine, the message INSERT PRO- 
GRAM TAPE, PRESS RETURN will be displayed. 


2. When the program tape is properly inserted, the message IS DATA STORED ON 
PROGRAM TAPE will appear. Type YES or NO and press RETURN. If your reply 
was YES, steps 3 and 4 below do not apply. 

3. Because the input data file is not on the program tape, the message INSERT DATA 
TAPE, PRESS RETURN will be displayed. 

4. When the data tape has been inserted, HEADER FILE £ will appear on the screen. 
Enter the Header File number and press RETURN. 

5. Data from the tape will be automatically entered and processed until the ALL DATA 
ENTERED message is displayed. If additional data is to be entered, follow the in- 
structions describing the ADD DATA key function. Otherwise, any of the function 
keys may be used to analyze the data. 

ADD DATA 


ADD DATA (key No. 2) allows data to be added to previously entered files. If the original 
entry was from the keyboard and was saved on tape, the added data will also be saved on 


tape. 


Ts 


If the original data came from a tape file, the added data will always be appended. 


Press ADD DATA. If the alarm bell sounds, a message on the display screen will in- 
dicate that a magnetic tape must be inserted. Insert the correct tape and press 
RETURN. 


. ADD DATA TO SAMPLE will be displayed. Enter X or Y and press RETURN. 


. To add data to sample X, follow the instructions beginning with step 5 of the ENTER 


DATA key description. For sample Y, begin with step 6. 


. When the added data points have been entered, they will be included in all further 


analyses. 


DELETE POINT 


DELETE POINT (key No. 3) is used to remove single observations from either sample. The 
deletion will eliminate the observation from internal storage, the tape file, and calculations. 


1. 


Press DELETE POINT (key No. 3). If the alarm bell sounds, insert the magnetic tape 
indicated by the message on the screen; press RETURN. 
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2. When the correct tape has been inserted, the following message will be displayed: 


FROM SAMPLE. Enter X or Y; press RETURN. Enter the value to be deleted; press 
RETURN. 


. The sample is searched and the observation is deleted from files and calculations. The 


message CORRECTION DONE indicates completion of the operation. If the data 
point is not within the designated sample, DATA NOT FOUND is displayed and no 
deletion occurs. 


CHANGE POINT 


CHANGE POINT (key No. 4) permits an incorrectly entered observation to be changed in 
all calculations. If the data is in storage, the stored data will also be corrected. 


1. Press CHANGE POINT. If the alarm bell sounds, go to step 2 below. If the bell does 


not sound, go directly to step 3. 


. Insert the correct tape indicated by the message on the screen; press RETURN. If the 


alarm bell sounds again, repeat this procedure. 


. The message IN SAMPLE is displayed. Enter X or Y, and press RETURN. Enter the 


incorrect value and the correct value, using a comma as a delimiter between the two 
values. Press RETURN. 


. The sample specified is searched for the incorrect value, and corrections are made in 


storage and in calculations. CORRECTION DONE is displayed to indicate completion. 
DATA NOT FOUND indicates that the value designated as incorrect was not found. 


. When this operation has been completed, other function keys may be used for analysis 


of the corrected data. 


LIST DATA 

When LIST DATA (key No. 5) is pressed, both data samples are listed on the screen. If a 
tape insertion is required, the alarm bell and the display screen will indicate the correct 
tape. If all data has not been stored, the message DATA IS NOT AVAILABLE TO BE 
LISTED will appear. 
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TWO SAMPLE t-TEST 


HISTO 


1. Press HISTO (key No. 6).If the alarm sounds, the message displayed will indicate a 
tape insertion requirement. Insert the tape and press RETURN. If the alarm sounds 
again, repeat this procedure. 


2. If the input data has not been stored, the following message will appear: HISTOGRAM 
MAY NOT BE CALCULATED. Otherwise, the message HISTOGRAM FOR SAMPLE 
will be displayed. Enter X or Y and press RETURN. 


3. The data minimum, maximum, and range for the selected sample are displayed with 
the following message: 


SELECT OPTION 
1 SPECIFY FOR CELLS 
2 SPECI FY CELL WIDTH 
3 SPECIFY OFFSET AND CELL WIDTH 


ENTER OPTION 


4. After entering the option number, press RETURN. The screen will then display either 
# OF CELLS; CELL WIDTH; or OFFSET, CELL WIDTH. Enter the appropriate in- 
formation and press RETURN. Remember, the number of cells for option 1 must not 
exceed 20. 


5. DO YOU WANT NORMAL CURVE? will now be displayed. If YES is entered, the 
curve will appear with the histogram when RETURN is pressed. For a NO answer, the 
histogram is drawn without the curve. 


pum | 


STATISTICS 


Descriptive statistics for both samples are displayed when STATISTICS (key No. 7) is pressed. 
Any required tape insertion is indicated by the alarm bell and a message on the display screen. 


Statistics displayed are: 
number of observations; 
sample mean, variance, and standard deviation; 
standard error of the mean; 
data minimum, maximum, and range; 
coefficient of variation; 
skewness; and 
kurtosis. 
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TWO SAMPLE t-TEST 


t-TEST AND CONFIDENCE INTERVALS 


t-TEST CI (key No. 8) causes calculation and printing of the t-test and confidence intervals 
for the means. If the alarm bell sounds when the key is pressed, a message on the display 
screen will indicate that a tape insertion is required. When the RETURN key is pressed, the 
t-test of equality of the means is calculated and printed with the confidence intervals for the 
means. 


C O C 


C 


STOP ENTRY 

The STOP ENTRY (key No. 12) is used only if data entry is not completed or if no imme- 
diate analysis is desired. Use of this key enables keyboard entered data to be stored for 
later additions or analysis. If the tape data storage option was selected, STOP ENTRY will 
initiate storage to the tape file. 


C ES O 


C 


EXAMPLE 1 

In this example, taken from Reference 2, (pages 49-53), IQ scores are obtained from 10 
girls whose parents have a high income and from 9 girls whose parents have a low income. 
The hypothesis states that the two sets of girls have the same mean IQ. 


Sample X is high income IQ measurements of 124, 114, 115, 106, 84, 96, 106, 126, 124, 
and 116. Sample Y, low income IQ data is 113, 97, 108, 95, 105, 69, 113, 98, and 118. 


C Es CGC 


Because of the small sample size, data will not be stored on tape. The program is loaded, 
and data is entered as shown below: 


I 


TWO SAMPLE T-TEST 
DO YOU WISH TO STORE DATA ON TAPE? NO 


C= “AS 4: 


c 


C 
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Press LIST DATA (key No. 5). 


124.0000 
96.0000 


113.0000 
69.0000 
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DATA 


X SAMPLE 
114.0000 
186.0000 


Y SRMPLE 
97.0008 
113. 60608 


ENTER DATA 


115.0008 
126.0000 


108.0000 
38.0680 
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106.0080 
124.0000 


95.0008 
118, 0090 


TWO SAMPLE t-TEST 


94.0820 
116.0006 


185.8008 
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TWO SAMPLE t-TEST 


Press HISTO (key No. 6) to display a histogram for both samples. For the X sample, 
specify 10 cells with a normal curve overlay. 


HISTOGRAM 
HISTOGRAM FOR SAMPLE ? X 


DATA MIN = 84 
DATA MAX = 126 
DATA RANGE = 42 


SELECT OPTION 
1 SPECIFY % OF CELLS 
2 SPECIFY CELL WIDTH 
3 SPECIFY OFFSET AND CELL WIDTH 


ENTER OPTION! 
& OF CELLS 10 


DO YOU WANT NORMAL CURVE YES 


C QZ G oL o C To 


N % REL FREQ X SAMPLE 
$ 58.0 


ae ek be ee ae ae DE 


C 


( it 


+8.840 +0.924 *1.088 *1.892 +1.176 
11012 
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For the Y sample, again specify 10 cells with a normal curve overlay. 


HISTOGRAM 
HISTOGRAM FOR SAMPLE ? v 


DATA MIN = 69 
DATA MAX = 118 
DATA RANGE = 49 


. SELECT OPTION 
1 SPECIFY 4 OF CELLS 
2 SPECIFY CELL WIDTH 
3 SPECIFY OFFSET AND CELL WIDTH 


ENTER OPTION1 
& OF CELLS 19 


DO YOU WANT NORMAL CURVE YES 


PN N x REL FREQ Y SAMPLE 
| 55.6 


due ren | 
a 


) 


+0.690 +0.788 +0.886 *8.984 *1.082 
11812 
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N= 

MEAN = 
VARIANCE = 
STD DEV = 
DATA MIN = 
DATA MAX = 
DATA RANGE = 


STANDARD 
ERR OF MEAN = 


COEFFICIENT 
OF VARIATION = 


SKEWNESS = 
KURTOSIS = 


T TEST OF EQUALITY OF MEANS = 


DESCRIPTIVE STATISTICS 


X SAMPLE Y SAMPLE 
18 9 
111.1 181.777777778 
179.655555556 215.194444445 
13. 4035650316 14, 6695073607 
84 69 
126 118 
42 49 
3.3331666625 3, 36283243343 
12, 0644149699 14.4132713653 
0. 768872598457 -1.1894220069 
2.6708246708234  1748,.44015146 


POOLED VARIANCE =196.379738562 


T-TEST 


WITH 17 DEGREES OF FREEDOM 


SIGNIFICANCE LEVELS : 
ONE TAIL TEST =0.08936 


X MEAN =111.1 


1.4478233179 


TKO TAIL TEST =0.17872 
CONFIDENCE INTERVALS FOR MEANS : 


LEVEL LOWER LIMIT 
% (102.328251373 , 119.871748627> 


967 


UPPER LIMIT 


95%  (181.388112173 , 120,€90887827> 
995 (97,322858778 + 124.877141222> 


Y MEAN 2181.777777? 


3-30 ` 


?8 
LEUEL LOMER LIMIT 


UPPER LIMIT 


904 (92.6825762195 , 119,2872979336> 
95% (98,4987615124 , 113.956794043> 
99% (85.3679996488 , 115.137595947) 


WMERN - YMEAN =9.32222222222 
LEV 


EL LOKER LIMIT 


UPPER LIMIT 


96% (-1.88137014212 . 28.5259145866) 
95% (-4.26576772416 , 22.9192121686> 
99% (-9.34239386871 , 27.9868383132) 
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TWO SAMPLE t-TEST 


The t-test for equality of means is 1.45 with significance level .18. Because of the size of 
this significance level, we do not reject the hypothesis of equality of the means. 


An inspection of the confidence intervals for the difference of the means shows that the 
confidence intervals contain both positive and negative values, (i.e., zero). The zero value 
again leads to the conclusion that we cannot reject the hypothesis of equal means. 


Although the difference of means, 9.32, is approximately 10% of the sample means, we do 
not reject equality of the means because of the large variance within the sample. The vari- 
ance is graphically illustrated by the histograms. 
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PAIRED t-TEST 


DESCRIPTION 


This program performs the paired t-test of population means equality for paired sample 
data. Data pairs may be entered from the keyboard or from magnetic tape files, and data 
may be corrected, deleted, or updated. 


Optional output may be selected from the following: 


. list of data pairs and differences; 
. histogram of either sample data, or of the differences; 
. data statistics including variance, standard deviation, skewness, and kurtosis; 


. t-test of equality of the means; 


an FP UU N = 


. confidence intervals for either sample or for the difference of the means. 


HARDWARE REQUIREMENTS 


A Tektronix 4050-Series Graphic System is required. If the Data Communications Interface 
(Option 1) is installed, minimum acceptable memory configuration is 16K bytes. 


[— ES . E E LL Es SEI C. 


C 


The number of data pairs which may be stored internally depends upon the available 
memory size of your machine. The following table shows the maximum number of data 
pairs permitted in internal storage without program modification . 


MEMORY SIZE MAXIMUM NUMBER OF 
DATA PAIRS 


E CES E SE 


For a 32K machine, up to 800 data pairs may be stored internally by changing the value of 
variable WO to 800. This variable is located on line 110 of program file 34. If the value of 
the variable is changed, data files on tape must then be updated to accomodate the change. 
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TAPE STRUCTURE 


The two types of taped files recognized and created by this program are Header Files 
and Data Files. Both must be in binary format. 


Header Files consist of two numbers; the Header File designation number, and the 
number of Data Files in the data set. 


Data Files each contain two types of information. The number of data observations 


contained in the file is followed by the data pairs in the form X4 Y4,X59Yo,. Xn Yn: 


The paired t-test file structure is illustrated below: 


HEADER FILE DATA FILE DATA FILE 


: Number of Number Data Data 
First Data Files in of Pair Pair 
File Number Data Set Pairs X4 Yı | X2 Yo Xn Yn 


METHODS 
The difference (D) of any XY data pair is X-Y. Thus, Dy = X, —Y,; Da = X, — Yp ete 


The paired t-test is the one sample t-test applied to the differences. Standard formulae used 
to calculate the descriptive statistics are shown in the Appendix. 


REFERENCES 
1. Dixon, W. J., Massey, F. J., /ntroduction To Statistical Analysis, New York: McGraw- 
Hill, 1957. 


2. Dunn, O. J., Clark, V. A., Applied Statistics: Analysis of Variance and Regression, 
New York: John Wiley and Sons, 1974. 


3. Peizer, D. B., Pratt, J. W., "A Normal Approximation for the Binomial, F, Beta, and 
Other Common Related Tail Probabilities 1“, Journal of the American Statistical 


Association, 63, 1968. 


4. Snedecor, G. W., Cochran, W. G., Statistical Methods, Sixth Edition, Ames lowa: lowa 
State University Press, 1967. 
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VARIABLES i 
VARIABLE 

A Data matrix [| 
Al Scratch 
A2 Scratch 
B1 Scratch L 
B2 Scratch 
B3 Scratch 
C Save summary statistics i 
D Cell width for histograms 
D1 Scratch 
D2 Scratch U 
D3 Scratch 
D4 Scratch j 
D5 Scratch U 
D6 Scratch 
D7 Scratch / L 
E$ String variable for alpha information V 
E Matrix for storing histogram data 
E1 Flag for saving data on tape 


O = do not save 


1 = save on data tape 

—1 = save on separate tape 
E2 Flag for tape type 

O = data tape in 

1 = program tape in 
E3 File required for overlaying program 
E4 Flag for step entry 


O = not performed 
1 = step entry performed 


E [— E TS ES ER 


E5 Program file currently in memory 

E8 Overlay key pressed 

F Beginning file number, number of files used 
| Scratch 

I1 Scratch 

J Scratch 

K Scratch 
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| PAIRED t-TEST 
VARIABLE USAGE 
n N Number of observations in the current file 
P Scratch 
Q Scratch 
I S Scratch 
T Scratch 
T1 Scratch 
i W2 Scratch 
WO Dimension of the data matrix A. 
X Scratch 
i X1 Scratch 
S2 Scratch 
X3 Scratch 
1 X4 Scratch 
Y1 Scratch 
n A Y2 Scratch 
Y3 Scratch 
Y4 Scratch 
| Plot 50: Stat Vol 1 3 -35 
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FUNCTION KEYS 


T-TESTS 


ES TIRE SERT 
SHIFT KEYS 


DATA STOP 
FROM TAPE ENTRY 


DELETE 
POINT 
She AUSSEN Ue ER 


SHIFT KEYS 


— PN334-2691-00 


HISTO i STATISTICS 


L4 


KEY KEY FUNCTION 

1 ENTER DATA allows entry of data from keyboard. 

2 ADD DATA causes data to be added to an existing file. 

3 DELETE POINT permits data pair to be deleted from memory and tape. 

4 CHANGE POINT allows data value to be changed. 

5 LIST DATA lists data previously entered and data pair differences. 

6 HISTO causes histogram of data to be displayed. 

7 STATISTICS calculates and displays the number of observations, the 
sample mean, sample variance, sample standard deviation, 
standard error of the mean, data minimum, data maxi- 
mum, data range, coefficient of variation, skewness, and 
kurtosis. 

8 T-TEST CI calculates t-test and confidence intervals for the mean of 
each sample and for the mean of the differences. 

11 DATA FROM TAPE | allows entry of data from taped file. 

12 STOP ENTRY terminates data entry from keyboard and permits storage 


of data, if tape storage option was elected. 
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PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the 
keyboard. 


AUTOMATIC LOADING FROM THE DIRECTORY 


1. 
2. 


Insert the program tape. 


Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


. Enter 25, press RETURN. 


4. When the I/O light goes out, execute the program. 


MANUAL LOADING FROM THE KEYBOARD 


. Insert the program tape. 


Enter FIND 26, press RETURN. 
Enter OLD; press RETURN. 
When the I/O light goes out, execute the program. 


PROGRAM EXECUTION 

If input data is to be entered from an existing magnetic tape file, follow the instructions 
concerning the DATA FROM TAPE (key No. 11). Otherwise, data must be entered from 
the keyboard as described below. 


ENTER DATA 


1. 


Plot 50: Stat Vol 1 


Press ENTER DATA (key No. 1). The screen displays DO YOU WISH TO STORE 
DATA ON TAPE? Reply YES or NO through the keyboard and press RETURN. If 
you replied NO, omit steps 2, 3, and 4; go to step number 5 below. If you replied 
YES, continue the instruction sequence. 


. DO YOU WISH TO STORE DATA ON PROGRAM TAPE? is displayed. If the file 


is on the program tape, type YES; press RETURN, and go to step 4. If the data is on 
a separate data tape, reply NO and press RETURN. 


. To establish the identity of the input data file, the message HEADER FILE is dis- 


played. Enter the correct Header File number and press RETURN. 
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4. If the alarm bell sounds, insert the tape indicated by the message on the screen and 
press RETURN. 


5. The display screen shows: 
ENTER DATA 
X,Y 


Enter the X and Y value, and press RETURN. Repeat this step for each pair. When 
all value pairs have been entered, use the function keys to perform data analysis. 


DATA FROM TAPE 
1. To allow entry of data pairs from an existing magnetic tape file, press DATA FROM 


TAPE (key No. 11). 


Data tapes used by this program must conform to the tape standards explained in 
the DESCRIPTION section. 


If the program tape has been removed from the machine, the message INSERT PRO- 
GRAM TAPE, PRESS RETURN will be displayed. 


2. When the program tape is properly inserted, the message IS DATA STORED ON 
PROGRAM TAPE will appear. Type YES or NO and press RETURN. If your reply 
was YES, steps 3 and 4 below do not apply. 


3. Because the input data file is not on the program tape, the message INSERT DATA 
TAPE, PRESS RETURN will be displayed. 


4. When the data tape has been inserted, HEADER FILE will appear on the screen. Enter 
the Header File number and press RETURN. 


5. Data from the taped file will be automatically entered and processed until the ALL 
DATA ENTERED message is displayed. If additional data is to be entered, follow 
the instructions describing the ADD DATA key function. Otherwise, any of the 
function keys may be used to analyze the data. 
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ADD DATA 


ADD DATA (key No. 2) allows data to be added to previously entered files. If the original 
entry was from the keyboard and was saved on tape, the added data will also be saved on 
tape. If the original data came from a tape file, the added data will always be appended. 


1. Press ADD DATA. If the alarm bell sounds, a message on the display screen will in- 
dicate that a magnetic tape must be inserted. Insert the correct tape and press 
RETURN. 


2. ADD DATA X,Y will be displayed. Enter the X,Y pair and press RETURN. Repeat 
this step for each pair to be added. 


3. When all pairs have been added, use the function keys to perform data analysis. 


DELETE POINT 


DELETE POINT (key No. 3) is used to remove data pairs from internal storage, calcu- 
lations, and from tape files. 


1. Press DELETE POINT. If the alarm bell sounds, insert the magnetic tape indicated by 
the message on the screen. Press RETURN. 


2. When the correct tape has been inserted, the following message will be displayed: 
DELETE X,Y. Enter the X,Y pair to be deleted and press RETURN. 


3. The data pair is deleted from files and calculations. 1f the message DATA NOT FOUND 
is displayed, no deletion has occurred. 


CHANGE POINT 


CHANGE POINT (key No. 4) permits an incorrectly entered pair to be changed in al! calcu- 
lations, storage, and tape files. 


1. Press CHANGE POINT. If the alarm sounds, go to step 2 below. If not, omit the next 
step and go to step 3. 


2. Insert the correct tape indicated by the message on the display screen. Press RETURN. 
Repeat this procedure if the alarm bell rings again. 


3-39 


Plot 50: Stat Vol 1 


TEST AND REGRESSIONS 


PAIRED t-TEST 


3. CHANGE X,Y TO X,Y is displayed. Enter the incorrect pair followed by the correct 
pair. Commas may be used as delimiters between the 4 values entered. Press RETURN. 


4. The data is searched for the incorrect value, and corrections are made in storage and 
in calculations. CORRECTION DONE is displayed to indicate completion. DATA NOT 


FOUND indicates that the value designated as incorrect was not found. 


5. When this operation has been completed, other function keys may be used for analysis 


of the corrected data. 


LIST DATA 


To list data previously stored internally or on magnetic tape, press LIST DATA (key No. 5). 
If a tape insertion is required, a reminder message will be displayed and the alarm bell will 


sound. Insert the correct tape and press RETURN. 


If the data has not been stored, DATA NOT AVAILABLE TO BE LISTED will be displayed. 


Remember, if the screen becomes full, no more data will be listed until HOME PAGE is 


pressed. 


HISTO 


HISTO (key No. 6) generates and displays histograms of X or Y data points, or of the 
differences. Necessary tape insertions are indicated by the alarm bell and a displayed 


message. The RETURN key must be pressed after the requested tape has been inserted. 


I f the data has not been stored, HISTOGRAM MAY NOT BE CALCULATED is displayed. 
Otherwise, the message HISTOGRAM FOR SAMPLE? is displayed. Enter either X, Y, or 


D (as in difference) to indicate the data to be used. 


The data minimum, maximum, and range will be displayed with the following message 


SELECT OPTION 
1 SPECIFY FOF CELLS 
2 SPECIFY CELL WI DTH 
3 SPECIFY OFFSET AND 


CELL WIDTH 
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After entering the option number, press RETURN. The next message will be # OF 
CELLS: or CELL WIDTH; or OFFSET, CELL WIDTH. Enter the information requested 
and press RETURN. If option 1 was selected, the number of cells specified must not ex- 
ceed 20. 


DO YOU WANT NORMAL CURVE will now be displayed. If your reply is YES, a normal 
curve overlay will appear in addition to the histogram. 


STATISTICS 


To calculate and display the descriptive statistics for each sample and for the difference, 
press STATISTICS (key No. 7). 


If the alarm bell sounds, insert the tape indicated by the message on the screen and press 
RETURN. 


t-TEST and CONFIDENCE INTERVALS 


t-TEST CI (key No. 8) initiates calculation and display of the t-test of the equality of means, 
and of the confidence intervals. Tape insertion requirements are indicated by the alarm bell 
and a message on the display screen. 


STOP ENTRY 


STOP ENTRY (key No. 12) is only used to store data for later analysis or for further entries 
at some other time. If all data has not been entered or if analysis is not wanted until later, 
press the STOP ENTRY key to initiate data storage. 


EXAMPLES 


Two new types of rations are fed to pigs. Which one will produce a higher weight gain? 


Two litter mates are selected from 12 litters and one member of each pair is randomly as- 
signed to receive ration A. The other will receive ration B. Weight gains are shown below: 


RATION A 30 31 34 29 26 32 31 35 32 38 34 29 
RATION B 31 26 24 28 29 30 28 29 32 32 26 29 
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Load the program and press ENTER DATA (key No. 1). Because the number of observations 
is small (12), we choose to store our data in memory rather than on tape. 


PAIRED SAMPLE T-TEST 
DO YOU MISH TO STORE DATA ON TAPE? NO 


c 
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LIST DATA (key No. 5) is pressed to display the data and differences. 


ENTER DATA 


X s Y 238,31 
DIFFERENCE =-1 
Ky Y 231,26 
led be 25 
f) 2343 
DIFFERENCE 210 
X , Y 229,28 
DIFFERENCE 21 
X , Y 226,29 
DIFFERENCE 2-3 
Y =32,38 
*OTFFERENCE 22 
X , Y 231,28 
DIFFERENCE =3 
Xs Y 235,29 
DIFFERENCE z6 
X , Y 232,32 
DIFFERENCE 29 
4, Y =38,32 


DIFFERENCE 36 
; Y 234,26 
DIFFERENCE 28 
K , Y 229,29 
DIFFERENCE =9 
X,Y 


DATA 

X Y DIFFERENCES 
38. 8808 31.8888 -1.0880 
31.0000 26.0000 5.0000 
34. 0880 24.0000 lo. 0009 
29.0000 28.0000 1.0000 
26.0000 29.0000 -3.0000 
32.0098 38. 8888 2.0000 

1.0000 28.0000 2.0890 
35.0000 29.0000 6.0800 
32. 8980 32.0090 0.0008 
38.0000 32.0008 5.0000 
34.0000 26.0008 8.0006 
29. 8008 29.0060 9.4080 
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PAIRED t-TEST 


For a histogram of the A ration data, press HISTO (key No. 6) and enter X. The data range 
is 12 so a cell width of 2 will provide 6 cells. Option 2 allows cell width specification. 


HISTOGRAM 
HISTOGRAM FOR SAMPLE ? x 
DATA MIN = 26 


DATA MAX = 38 
DATA RANGE = 12 


SELECT OPTION 
1 SPECIFY # OF CELLS 
2 SPECIFY CELL WIDTH 
3 SPECIFY OFFSET AND CELL WIDTH 


ENTER OPTION2 
CELL WIDTH 2 


DO YOU WANT NORMAL CURVE YES 


N x REL FREQ X SAMPLE 
9$. 41.7? 


> 


42.608 +2.800 42.000 +3.200 +3.400 +3.600 +3.800 


*1011 
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PAIRED t-TEST. 


The same procedure is followed for the Y sample (ration B). 
HISTOGRAM 
HISTOGRAM FOR SAMPLE ? Y 


DATA MIN = 24 
DATA MAX = 32 
DATA RANGE = 8 


SELECT OPTION 
1 SPECIFY % OF CELLS 
2 SPECIFY CELL WIDTH 
3 SPECIFY OFFSET AND CELL WIDTH 


ENTER OPTIONZ 
WIDTH 2 


i. E. C C C Ts 


CELL 
DO YOU WANT NORMAL CURVE YES 


- 
— 


N x REL FREQ Y SAMPLE 


? $8. 3 { 


p. 1. 


*2.400 +2.608 +2.800 +3.888 +3.200 
11011 


c 


EO UE. TS I T 
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To display a histogram of the differences, press HISTO and enter D. In this example, 
option 3 is chosen with an offset of —4, and cell width of 3. 


HISTOGRAM 
HISTOGRAM FOR SAMPLE ? D 


DATA MIN = -3 
DATA MAX = 10 
DATA RANGE = 13 


SELECT OPTION 
1 SPECIFY 4 OF CELLS 
2 SPECIFY CELL MIDTH 
3 SPECIFY OFFSET AND CELL WIDTH 


ENTER OPTION3 
OFFSET, CELL WIDTH -4,3 


DO YOU WANT NORMAL CURVE YES 


N x REL FREQ DIFFERENCES 
6, 58.8 


~0. 488 ~@. 188 +8. 208 +9.508 +8, 888 
31011 
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The STATISTICS (key No. 7) will initiate the display shown below. 


DESCRIPTIVE STATISTICS 


X SAMPLE Y SAMPLE 
N= 12 

MEAN = 31.75 28. 6666666667 
VARIANCE = 18. 2045454545 6.06060606061 
STD DEV = 3. 19445542378 2.46182981959 
DATA MIN = 26 24 

DATA MAX = 38 32 

DATA RANGE = 12 8 

ERR GF MEAN = 0. 922159949418 0. 710663054313 
OF VARIATION = 10,0612769253 8.587 77844042 
SKEWNESS = 8. 184598883872 ~8, 323854985825 
KURTOSIS = 2. 7678285326 2. 31090906007 

T-TEST 


DIFFERENCES 


3.08333333333 
15.5378787879 
3.94181161243 
-3 
18 
13 


1.13798299777 


127.842538781 


6.223651408036 
2. 08433453988 


T-TEST OF MEAN DIFFERENCES = O IS 2.70966272115 


WITH 11 DEGREES OF FREEDOM 


SIGNIFICANCE LEVELS : 
ONE TAIL TEST 28.0873? 
TWO TAIL TEST =0.01474 
MEAN DIFFERENCE = 3.68333233333 


CONFIDENCE INTERVALS FOR MEAN DIFFERENCE : 


LEVEL LOWER LIMIT UPPER LIMIT 


90x (1. 03926949044 s 5.127 


37815622) 


95% (0.578174982564 , 78 S884916641) 
99%  (-8.451459486936 , 6.6181261536) 


The t-test of equality of means is 2.71, with a two-tail significance level .015. The null 


hypothesis is rejected at the .05 level. The 95% confidence interval for the mean differences 


3-46 i contains only positive numbers, and the null hypothesis is also rejected at the 596 level. 
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DESCRIPTION 


This program will fit any of the following functions to data of the form (X,Y): Y = AX; 
Y = A+BX; Y = AeB X; Y = 1/(A+BX); Y = A+B/X; Y = A+B log X; Y = AXE; and 
Y = X/(A*BX). 


Data, entered from the keyboard or from magnetic tape, may be added, deleted, corrected, 
or listed. 


Data can be plotted, the best fitting curve selected on the basis of minimum residual error, 
maximum R-square, or smallest maximum absolute deviation, and a specified range of data 
values can be used to fit curves or plot data. 


Each curve can be plotted against data, and the regression coefficients A and B can be calcu- 
lated and displayed. Maximum absolute deviation as well as Y and deviation or residual 
for al! values of X can be calculated. 


For the curve Y = A*BX, 95% confidence bands for prediction of Y may be plotted. 


The program will also calculate the R-square, and will estimate Y for any value of X entered 
from the keyboard. 


HARDWARE REQUIREMENTS 


A Tektronix 4050-Series Graphic System with a minimum 16K Read/Write memory is 
required. 


All data must be stored internally or on tape when the data is to be listed, plotted, trimmed, 
plotted against a curve, or used to calculate goodness of fit. The number of data pairs which 
may be stored internally is dependent upon machine memory size as shown by the table 
below: 


MEMORY SIZE MAXIMUM DATA PAIR CAPACITY 


16K 75 pairs 
24K 75 pairs 
32K 300 pairs 


Data sets larger than allowable for internal storage must be stored on magnetic tape files. 
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TAPE STRUCTURE U 
There are two versions of this program on the program tape. The first version, located in 
files 40, 41, 42, will run on a machine with 16K memory, but wil! require program over- i 
lay. Version two is totally contained in file 43 and may be executed on 24K or 32K | 
machines without program overlay. 

å 
When file 40 is loaded, directly or indirectly, and either the ENTER DATA key or DATA L 
FROM TAPE key is pressed, the program will be fully loaded, automatically. 
All data tapes used by this program must be binary coded, and each must contain a Header | 
File and one or more data files. The Header File, used to load data into memory, consists of 
the file number assigned to the first data file, and the number of files in the data set. T 


3-48 


The initial entry in each data file is the number of data pairs in the file. This is followed 
by the data pairs in the form Xi Yq: X4, Y», . .. When data is stored on data tapes, mes- 1 
sages on the display screen will prompt insertion of the correct tape. 


When large amounts of data are required, the data file should be stored on magnetic tape | Y 
for efficient program execution. 


C 


METHODS 


Each of the non-linear equations is transformed to an equation which is linear in parameters 
A and B. The least squares fit is calculated to the transformed data, and the parameters are 
then back-transformed to the original equation. 


The transformations are: 


ORIGINAL EQUATION TRANSFORMED EQUATION RESTRICTION 


Y=AX unchanged none 

Y = AtBX unchanged none 

Y = AeBX In Y = In A+BX Y>0 
Y =1/(A+BX) 1/Y = A+BX Y>0 
Y=A+B/X unchanged X>0 

Y = A+B log X unchanged X20 
Y = AXB In Y = In A+B In X XY >0 
Y = x/(A+BX) 1/Y= A/X*B X,Y >0 


Goodness of fit measures used by this program are: 


C 
Se ecc E 


1. R-square, an index determination; 
2. residual error, often called standard error of the estimate; and 


3. maximum absolute residual. 
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To define these measures, let Y = g(X) be any of the 8 equations; let (X; Y;) be the data. 


R-square is defined as: 


Z(Y; —9(X,))? 
R-square |— ————————- where Y = ZY./N 
Z(Y; YE 


Z(Y, —g(X;))? 
Residual error is: ————————— 
N —2 


Maximum absolute residual is defined as max (IY; — g(X)). 


With this definition, R-square may be negative. A negative R-square for a particular equation 
indicates that the equation is not appropriate for the data. 


When the SELECT BEST FIT key is pressed, and when all data is available, R-square, re- 
sidual error, and maximum absolute deviation are calculated and displayed for each equation 
that may be fit. A message indicates which equation has the smallest maximum absolute 
residual, and which equation has the maximum R-square. 


2 


REFERENCE 
Draper, N., Smith, H., Applied Regression Analysis, New York: John Wiley and Sons, 1966. 


VARIABLES 


VARIABLE USAGE 


Data array 
A1 Scratch 
A2 Scratch 
A3 Scratch 
A4 Scratch 
A5 Scratch 
B Coefficient matrix 
B1 Flag = 1 if coefficient calculated 
B2 Flag = —1 if equations 3, 4, 7, 8 can not be fit 
B3 Flag = —1 if equations 5, 6, 7, 8 can not be fit 
B4 Flag = 4 to change range of plot 
a B5 Flag = 4 to change range of fit 


C Summary statistics matrix 
D Residual error and maximum residual matrix 
Plot 50: Stat Vol 1 D3 Scratch 3 -49 
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VARIABLES 


VARIABLE 
D4 
E1 


E2 
E3 
E4 


AS =o 


oozr 


= 
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USAGE 

Scratch 

Flag = —1 data on separate data tape; 0 data not on tape; 
1 data on program tape 

Flag = Odata tape in; 1 program tape in 

Scratch 

Flag = O stop entry not performed; 1 stop entry per- 
formed 

Overlay key pressed 

Flag = 1 all data in memory; O all data not in memory 

String variable for alpha information 

Header File number, number of files used 

Array for drawing axes 

Scratch 

Scratch 

Scratch 

Scratch 

Number of data pairs 

Scratch 

Scratch 

Equation number to plot, etc. 

Scratch 

Dimension of data array 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 
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SIMPLE REGRESSIONS 


DATA 
FROM TAPE 


ENTER DATA 


ADD DATA 
DELETE POINT 


CHANGE POINT 
LIST DATA 


1 
A+BX 


A+B/X 


A+BlogX 


AXB 
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Kr 


SHIFT KEYS 


STOP TRIM 
ENTRY DATA 


DELETE CHANGE 
POINT POINT 


SHIFT KEYS 


SELECT 
BEST FIT AX 


A+ BLOGX 


FUNCTION 


allows entry of data from keyboard to memory and causes 
data to be stored to tape when input was from tape file. 


allows additions to data previously entered. 

deletes a designated data point from storage, tape file, 
and calculations. 

permits correction of a data pair. 

lists data previously stored. 


produces a plot of regression curve against the data, es- 
timation of Y for X, and residual calculation for equation 
Y =1/(A+BX). 

produces a plot of regression curve against the data, es- 
timation of Y for X, and calculation of residuals for 
equation Y = A+B/X. 

produces a plot of regression curve against the data, es- 
timation of Y for X, and calculation of residuals for 
equation Y ~ A+BlogX. 

pre uces 5 tot of regression curves against the data, es- 
timatic:: of Y for X, and calculation of residuals for 
equation Y = AXB. 
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A+BX 
DATA FROM 
TAPE 
STOP ENTRY 


TRIM DATA 


PLOT DATA 


SELECT BEST 
FIT 


AX 


PROGRAM LOADING 


FUNCTION 


produces a plot of regression curves against the data, es- 
timation of Y for X, and calculation of residuals for 
equation Y = X/(A+BX). 


allows entry of data from data tape or program tape. 


terminates data entry and saves data for later analysis or 
addition of more data. 


allows execution to be restricted to a portion of data pre- 
viously entered. 


initiates plot of data from taped data file or from internal 
storage. 

lists coefficients for equations and descriptive statistics for 
data in memory or tape files; calculates R-square. 


permits selection of plot of regression curve against data; 


calculates residuals; estimates Y for X for equation Y = AX. 


produces a plot of regression curve against the data, 
estimates Y for X, calculates residuals for equations 
Y=A+ BX. 


produces a plot of regression curve against the data, 
estimates Y for X, calculates residuals for equation 
Y = AeBX, 


This program may be loaded automatically from the directory, or manually from the key- 
board. The STATISTICS, VOLUME 1 program tape is required. 


AUTOMATIC LOADING FROM THE DIRECTORY 


1. Insert the program tape. 
2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 


THE PROGRAM NUMBER YOU WANT. 


3. Enter 26; press RETURN. 


4. When the 1/0 light goes out, execute the program. 
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MANUAL LOADING FROM THE KEYBOARD 
1. Insert the program tape. 
2. Enter FIND 40; press RETURN. 
3. Enter OLD; press RETURN. 
4. When the I/O light goes out, execute the program. 


mo ee fo Jess —J 


PROGRAM EXECUTION 


Program execution procedures are determined by the input method and data storage needs. 
The 5 options from which one selection must be made are as follows: 


1. Input data from the program tape; 

2. Input data from a previously created data tape; 

3. Input from the keyboard, data to be stored on the program tape; 
4. Input from the keyboard, data to be stored on a data tape; and 


) 


5. Input from the keyboard, no data to be stored on tape. 


ud: LAE II. ER 


Execution procedures for each of these options are described on the following pages. 


OPTION 1, EXECUTION WITH DATA INPUT FROM THE PROGRAM TAPE 


Execute the program using a data set previously stored on the program tape. 


1. Press DATA FROM TAPE (key No. 11). IS DATA ON PROGRAM TAPE is displayed. 


2. Enter YES from the keyboard; press RETURN. 
When all data from the program tape has been read into memory, the I/O light 
goes out and ALL DATA ENTERED is displayed. 


3. Use the appropriate function keys to execute program operations. 


OPTION 2, EXECUTION WITH DATA INPUT FROM DATA TAPE 


To execute the program with input from an existing data tape: 


1. Press DATA FROM TAPE (key No. 11). IS DATA ON PROGRAM TAPE is displayed. 


n 2. Enter NO from the keyboard; press RETURN. 
HEADER FILE is displayed. 


q Plot 50: Stat Vol 1 
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3. Enter the numeric designation of the data file to be used. Press RETURN. The alarm 
bell sounds and INSERT DATA TAPE, PRESS RETURN is displayed. 


4. Insert the correct data tape and press RETURN. When the data transfer is complete, 
the 1/0 light goes out and ALL DATA ENTERED is displayed. 


5. Use the appropriate function keys to execute program operations. 


OPTION 3, EXECUTION WITH KEYBOARD INPUT, DATA TO BE STORED ON PRO- 
GRAM TAPE 


To execute the program and store keyboard entered data on the program tape: 


1. Press ENTER DATA (key No. 1). DO YOU WISH TO STORE DATA ON TAPE 
is displayed. 


ES Es. 


2. Enter YES from the keyboard; press RETURN. DO YOU WISH TO STORE DATA 
ON PROGRAM TAPE is displayed. 


C 


3. Enter YES from the keyboard and press RETURN. 


ENTER DATA 
X,Y 


is displayed. 
4. Enter the value of X and Y, separated by a comma. Press RETURN after each X,Y pair. 


5. When all data has been entered, use the appropriate function keys to perform other pro- 
gram operations. Data storage will be automatic. 


ES E ELLER: 


6. If no analysis is desired, press STOP ENTRY (key No. 12) to store the data and up- 
date the Header File. 


E 


OPTION 4, EXECUTION WITH KEYBOARD INPUT, DATA TO BE STORED ON A 
DATA TAPE 
To execute the program and store keyboard entered data on a data tape, follow the steps 
below. 
1. Press ENTER DATA (key No. 1). DO YOU WISH TO STORE DATA ON TAPE is 
displayed. 


=. 


C 


p. ES 
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. Enter YES from the keyboard and press RETURN. DO YOU WISH TO STORE DATA 


ON PROGRAM TAPE is displayed. 


. Enter NO from the keyboard and press RETURN.HEADER FILE is displayed. 


. Enter the number of the header file and press RETURN. INSERT DATA TAPE, 


PRESS RETURN is displayed and the alarm bell sounds. 


. Insert the correct data tape and press RETURN. 


ENTER DATA 
XY 


is displayed. 


. Enter the X, Y pair values with a comma as a value delimiter. Press RETURN after 


each pair. 


. When all data pairs have been entered, use the appropriate function keys to perform 


program operations. Data storage will be automatic. 


. If no data analysis is desired, press the STOP ENTRY key (No. 12) to store the data 


and update the Header File. 


OPTION 5, EXECUTION WITH KEYBOARD INPUT AND NO DATA STORAGE 
To execute the program without storing data for later use, proceed as follows: 
1. 


Press ENTER DATA (key No. 1). DO YOU WISH TO STORE DATA ON TAPE is 
displayed. 


Enter NO from the keyboard and press RETURN. 


ENTER DATA 
X,Y 


is displayed. 


. Enter the X and Y values, separated by commas. Press RETURN after each pair. 


4. When all data has been entered, use the function keys to perform any of the operations 
described on the next pages. 
1 
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PROGRAM OPERATIONS 


Upon completion of the initial data entries, operations are performed with the function 
keys. The operations described here may be performed in any sequence. 


ADD DATA 


ADD DATA (key No. 2) will permit addition of new data to the file currently in 
memory. If the file also resides on tape, the alarm bell sounds to indicate that the tape 
should be inserted. A message on the screen will identify the kind of tape required; pro- 
gram tape or data tape. Data entered will be added to the tape file automatically. 


When the tape has been inserted, press RETURN. 


ENTER DATA 
X,Y 


will be displayed, and values may be entered. The RETURN key must be pressed after each 
pair of values has been typed. 


Other operations may now be performed, or the STOP ENTRY (key No. 12) may be used 
to update the file without other operations. 


DELETE POINT 


DELETE POINT (key No. 3) initiates deletion of a data pair value. The values specified 
will be deleted from internal storage and from the tape file, when applicable. All subsequent 
calculations and displays will exclude the deleted pair. 


If the alarm bell sounds when DELETE POINT is pressed, insert the tape indicated by the 
message on the screen and press RETURN. 


DELETE X,Y will appear on the screen. Type the value pair to be deleted, and press 
RETURN. The message CORRECTION DONE indicates that the values have been deleted. 


ENTER DATA 
X,Y 


is displayed to allow other values to be entered. 


If the pair selected for deletion is not found, DATA NOT FOUND is displayed and all pro- 
cessing stops. 
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CHANGE POINT 


If an incorrect value has been entered, CHANGE POINT (key No. 4) provides a method 
of correcting the data. When the key is pressed, any necessary tape insertion is indicated by 
the alarm bell and a message on the screen. 


The screen displays the message CHANGE: INCORRECT X,Y TO CORRECT X,Y. Enter 

the incorrect values followed by the correct values, using a comma to separate each entry. 
When RETURN is pressed, CORRECTION DONE will indicate that the data has been changed 
internally and on tape, if applicable. Additional data may be added or other keys may be 

used to perform analysis. 


If the designated values are not found, all processing stops and DATA NOT FOUND is dis- 
played. Analysis may be resumed by pressing any other function key. 


LIST DATA 


Data currently in internal storage or residing on a tape file will be listed when the LIST 
DATA key (No. 5) is pressed. If a tape insertion is required, the alarm will sound and a 
message will appear on the screen. Listing of the data will continue until the entire data set 
has been displayed. The HOME PAGE key must be used to clear the screen when the display 
area is full. 


1 
A+BX 
This key causes calculation and display of the details for curve Y = 1/(A*BX). If the alarm 
bell sounds when the key is pressed, tape insertion requirements will be shown on the display 
screen. When the displayed instructions have been satisfied, the following message will ap- 
pear on the display screen: 


(Key No. 6) 


FOR EQUATION Y = A+BX SELECT OPTION 


1 PLOT EQUATION 
2 RESIDUALS 
3 ESTIMATE Y FOR X 


4 PLOT EQUATION AND CONFIDENCE BANDS 
ENTER OPTION 


Enter option number 1, 2, 3, or 4, and press RETURN. The appropriate operation will be 
performed and results will be displayed. If the operation is not completed, a displayed 
message will identify the reason. 


Ptot 50: Stat Vol 1 


3-57 


TEST AND REGRESSIONS 


SIMPLE REGRESSIONS . 


3-58. 


A+B/X (Key No. 7) 


This key causes calculation and display of the details for curve Y = A+B/X. If the alarm bell 
sounds when the key is pressed, tape insertion requirements will be shown on the display 
screen. When the displayed instructions have been satisfied, the following message will ap- 
pear on the display screen: 


FOR EQUATION Y = A*B/X SELECT OPTION 


1 PLOT EQUATION 

2 RESIDUALS 

3 ESTIMATE Y FOR X 
ENTER OPTION 


Enter option number 1, 2, or 3, and press RETURN. The appropriate operation will be per- 
formed and results will be displayed. If the operation is not completed, a displayed message 
will identify the reason. 


A+BlogX (Key No. 8) 

This key causes calculation and display of the details for curve Y = A*BlogX. If the alarm 
bell sounds when the key is pressed, tape insertion requirements will be shown on the dis- 
play screen. When the displayed instructions have been satisfied, the following message will 
appear on the display screen: 


FOR EQUATION Y = A+BlogX SELECT OPTION 


1 PLOT EQUATION 

2 RESIDUALS 

3 ESTIMATE Y FOR X 
ENTER OPTION 


Enter option number 1, 2, or 3, and press RETURN. The appropriate operation will be per- 
formed and results will be displayed. If the operation is not completed, a displayed message 
will identify the reason. 


AX? (Key No. 9) 


This key causes calculation and display of the details for curve Y = AXB. If the alarm bell 
sounds when the key is pressed, tape insertion requirements will be shown on the display 

screen. When the displayed instructions have been satisfied, the following message will ap- 
pear on the display screen: 
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FOR EQUATION Y = AXB SELECT OPTION 


1 PLOT EQUATION 

2 RESIDUALS 

3 ESTIMATE Y FOR X 
ENTER OPTION 


Enter option number 1, 2, or 3, and press RETURN. The appropriate operation will be per- 
formed and results will be displayed. If the operation is not completed, a displayed message 
will identify the reason. 


X 
AIBX (Key No. 10) 
This key causes calculation and display of the details for curve 
X 
Y = 
A+BX 


If the alarm bell sounds when the key is pressed, tape insertion requirements will be shown 
on the display screen. When the displayed instructions have been satisfied, the following mes- 
sage will appear on the display screen: 


X 
FOR EQUATION Y = SELECT OPTION 
A*BX 
1 PLOT EQUATION 
2 RESIDUALS 


3 ESTIMATE Y FOR X 
ENTER OPTION 


Enter option number 1, 2, or 3, and press RETURN. The appropriate operation will be per- 
formed and results will be displayed. If the operation is not completed, a displayed message 
will identify the reason. 


STOP ENTRY 

STOP ENTRY (key No. 12) should be used when keyboard entered data is to be stored 
to tape, and when no immediate analysis is desired. When the key is pressed, the last data 
file entered will be written to the specified tape file. If that file already exists on tape, the 
data will be added to it. The control information required for the Header File will also be 
written or updated. 


When the operation has been completed, REMOVE TAPE will be displayed. 3.59 
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TRIM DATA 


TRIM DATA (key No. 13) makes possible the change of plotting and listing range with- 
out altering the fit of the curves. This permits a close examination of the data and regression 
equations within the selected range, yet retains the original fit. 


The TRIM DATA may also be used to change the range of fit to a narrower range within 
the original data. 


When the key is pressed, the alarm and a message on the screen will indicate any required 
tape insertions. When the tape requirement has been satisfied, the message displayed will be: 


DO YOU WISH TO 


1. CHANGE PLOT RANGE ONLY 
2. CHANGE RANGE OF FIT AND PLOT 
3. RESET FIT AND PLOT RANGE TO ORIGINAL DATA 
4. RESET PLOT RANGE ONLY 
SELECT OPTION 


Enter the number of the option desired and press RETURN. 


I f option 1 was selected, CHANGE PLOT RANGE TO will appear on the screen. The 
minimum and maximum value of the range chosen should be typed. Å comma serves as å 
delimiter between the two values, and the RETURN key causes the values to be entered. 


The X data range, range of plot, and range of data fit are displayed, and plots, equations, 
and listings are based on them until further alterations are made. 


Option 2 causes CHANGE RANGE OF FIT TO to be displayed. When the minimum and 
maximum ranges have been entered through the keyboard, press RETURN. The summaries 
will be re-calculated to conform to the new data limits, and ALL DATA ENTERED will be 
displayed to indicate completion. Use other function keys to perform operations based upon 
the new fit. 


Options 3 and 4 are used to reset altered ranges to their original form. Data is re-entered 
from storage, and ALL DATA ENTERED is displayed to acknowledge the restoration. All 
subsequent operations and calculations will be based on the original data. 


If data is not available to permit the use of this key, OPERATION MAY NOT BE PER- 
FORMED will be displayed. 
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PLOT DATA 


PLOT DATA (key No. 16) provides a plot of available data. When the key is pressed, the 
plot will appear or tape insertion requirements will be signaled by the alarm bell and a dis- 
played message. 


SELECT BEST FIT 


To cause calculation of the coefficients for all functions that can be fit, press SELECT 
BEST FIT (key No. 17). Any necessary tape insertions will be indicated by the alarm bell 
and a displayed message. 


If data is available, R-square and standard error of the estimate will be calculated and dis- 
played. The maximum absolute residual is calculated for each regression that may be 
calculated. 


AX (Key No. 18) 


This key causes calculation and display of the details for curve Y = AX. If the alarm bell 
sounds when the key is pressed, tape insertion requirements will be shown on the display 
screen. When the displayed instructions have been satisfied, the following message will ap- 
pear on the display screen: 


FOR EQUATION Y = AX SELECT OPTION 


1 PLOT EQUATION 

2 RESIDUALS 

3 ESTIMATE Y FOR X 
ENTER OPTION 


Enter option number 1, 2, or 3, and press RETURN. The appropriate operation will be per- 
formed and results will be displayed. If the operation is not completed, a displayed message 
will identify the reason. 


A+BX (Key No. 19) 


This key causes calculation and display of the details for curve Y = A+BX. If the alarm bell 
sounds when the key is pressed, tape insertion requirements will be shown on the display 
screen. When the displayed instructions have been satisfied, the following message wil! appear 
on the display screen: 


FOR EQUATION Y = A+BX SELECT OPTION 


1 PLOT EQUATION 

2 RESIDUALS 

3 ESTIMATE Y FOR X 

4 PLOT EQUATION AND CONFIDENCE BANDS 
ENTER OPTION 
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Enter option number 1, 2, 3, or 4, and press RETURN. The appropriate operation will be 
performed and results will be displayed. If the operation is not completed, a displayed 
message will identify the reason. 


AeBX (Key No. 20) 

This key causes calculation and display of the details for curve Y = AeBX. |f the alarm bell 
sounds when the key is pressed, tape insertion requirements will be shown on the display 
screen. When the displayed instructions have been satisfied, the following message will ap- 
pear on the display screen: 


FOR EQUATION Y = AeBX SELECT OPTION 


1 PLOT EQUATION 

2 RESIDUALS 

3 ESTIMATE Y FOR X 
ENTER OPTION 


Enter option number 1, 2, or 3, and press RETURN. The appropriate operation will be per- 
formed, and results will be displayed. If the operation is not completed, a displayed message 
will identify the reason. 


EXAMPLE 
Find an equation which adequately represents the data below: 
X Y 
3 5.5 
4 7.75 
5 10.6 
7 13.4 
9 18.5 
12 23.6 
13 26.2 
14 27.8 
17 30.5 
20 33.5 
23 35 
25 35.4 
34 41 
38 42.1 
42 44,4 
45 46.2 
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SIMPLE REGRESSIONS 


Load the program and press ENTER DATA. Because the number of observations is small, 
data may be stored internally and there is no need to create a tape file. 


SIMPLE REGRESSIONS 


na 4. 3 


DO YOU WISH TO STORE DATA ON TAPE NO 


ENTER DATA 


wo ILE 3 


5 


45,46.2 


nn HM RR REN 
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List the data by pressing the LIST DATA key. 


SA meds Faca 
3 


X Y 

3 9.9 
4 7.75 
3 10.6 
? 13.4 
9 18.5 
12 23.6 
13 26.2 
14 27.8 
1? 38.5 
28 33.5 
23 35 
25 35.4 
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Press the PLOT DATA key to produce a plot of all untrimmed data. 1 


Ue Es 


There is a break point in the plot somewhere between X = 14 and X = 17. The data is nearly 
linear below 14 and above 17. 


C 
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SELECT BEST FIT (key No. 17) provides a table of criteria for all curves. 


SELECT BEST FIT 


EQUATION 
a B RES ERROR R-SQUARE MAY DEVIATION 
Y = Axx 
1.25433 56. 76697 9,70018 18.244936 
Y s A+ BiX 
3.95581 0.90726 20.24963 0.82385 7.17758 
Y = AXEXPCBXS) 
19,86252 9.03369 79.24698 8,528145 18.612680 
Y = A+ BEX) 
8.83887 -0,00227  14257.32734 -74.390152 341.45102 
Y = A+ BX 
48.61515 -133.32444 32.274898 8.82954 8.32633 
Y= A+ BtLOGCO 
| FN -14.18756 15. 61145 1.54151 8.99196 2.73892 
Y = AkXtB 
3.14259 8.75176 28.523394 8.39160 8. 76724 
= X/(R + BXX) 
8.49607? 8.08542 49.55325 8.74356 14. 60028 


EQUATION Y = A + BALOGCX) HAS MAXIMUM R-SQUARE 
EQUATION Y = À + BYLOGCX> HAS MINIMUM MAXIMUM ABSOLUTE RESIDUAL 


en en ee | 


The equation Y —14.18756 + 15.61145 log X has the highest R-square and the smallest 
maximum deviation. 


c uc 


The equation (Y = 1/(A*BX) has a negative R-square indicating that the equation form is 
not appropriate for this data. 


1 3 3 
2 
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Any of the individual function keys except Y = A+BX will provide the following display, 
using the equation corresponding to the key: 


FOR EQUATION Y = A*X SELECT OPTION 


1 PLOT EQUATION 

2 RESIDUALS 

3 ESTIMATE Y FOR X 
ENTER OPTION 


For the Y = A+BX key, the display will be: 


FOR EQUATION Y = A+BX SELECT OPTION 


1 PLOT EQUATION 

2 RESIDUALS 

3 ESTIMATE Y FOR X 

4 PLOT EQUATION AND CONFIDENCE BANDS 
ENTER OPTION 


All equations are plotted against the data, and for Y = A+BX, the 95% confidence bands 
are plotted. 
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Y = REN 


1.25433247394 


Å 


R-SQUARE = 
0. 708179499795 


3 3 


RES ERROR 
56. 7669678487 


MAXCABS CRESIDUAL > > 
10.2449613272 
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Y = ASEXPCBEX) 


Aa 
18. 8625246194 
Ba 
0.8396927121718 
R-SQUARE = 
0.581449014839 
RES ERROR 
79.2469815562 


MAXCABSCRESIDUAL >> 
18.6128821636 


Lo L— to te. 
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Y æ 1/(R + BEX) 


Å = 
0.0980729903147 


B = 
-0.088227336467726 


R-SQUARE = 
-74. 3015231463 
RES ERROR 
14257. 3273415 


MAX CABSCRESIDUAL > > 
341.451017642 
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Y = À + BX 


Å =z 
40.6151537502 
Bs 
7-133. 324443909 
R-SQUARE = 
8.829537338628 


RES ERROR 
32. 27488048264 
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Y * À + BLOGG) 


Az 
-14.1875642197 


15.6114538108 


R-SQUARE = 
8.991838335806 


RES ERROR 
1.54151439852 


MAXCABSCRESIDUAL > > 
2.7909221923? 
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Y = ARXIB 

Az 
3. 14258996212 

B= Ú 
0.751758295473 

R-SQUARE = Li 
0. 891600734877 


RES ERROR 
20.5239381869 


MAXCABSCRES IDUAL)) 
8.76724120826 | 
V . 
U 
U 
U 
U 
L 
i 
— 
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a oe, CF. 3 


Y = K/CR + BEX) 


Z3 


Å 
0.496065131274 


B 
0.00542362425122 


R-SQUARE = 
0.743561089159 


23 23 


RES ERROR 
48.5532475596 

MAXCABS (RESIDUAL >> 
14.6002757187 
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C 
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Y 2 f + BK 


As 
9.953580897533 


Bs 
0.907253741606 


G G t+ 


R-SQUARE = 
8. 893054784983 


RES ERROR 
20. 2486334192 


MAX CABSCRESIDUAL > > 
7.17758122351 


b= I E 


C 
CO (2 ba 


p Es 


The equation Y = A+B log X appears to be the best fitting equation. The extreme oscillation 
about X = 30 on the Y = 1/A+BX plot indicates that the equation is inappropriate. 


cC 


To illustrate the use of TRIM DATA (key No. 13) data is trimmed to the range X = 3 to 
X = 14, and the best fitting equation in this range is found. 


C 
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TRIM DATA 
RANGE OF X DATA IS 3 TO 45 


EQUATIONS HAVE BEEN FIT TO DATA WITH 
X RANGE FROM 3 TO 45 


THE PLOT RANGE IS 3 TO 45 


aol! E. 3. | 


DO YOU WISH TO : 
1 CHANGE PLOT RANGE ONLY 
2 CHANGE RANGE OF FIT AND F 
3 RESET FIT AND PLOT RANGE 10. ORIGINAL DATA 
4 RESET PLOT RANGE ONLY 


SELECT OPTION 


2 
CHANGE RANGE OF FIT TO 3.14 
ALL DATA EHTERED 


ced. uM a 


All data is automatically re-entered from storage or tape, and no further operation should 
be attempted until ALL DATA ENTERED is displayed. 


j 


If option 1 had been excercised (rather than option 2, as shown), equations in the 3 to 14 
range could have been closely observed without changing the fit. 


2J 23 


The LIST DATA key is used to verify that only data in the 3 to 14 range is being considered. 


OQ 


DATA 
n X Y 
' 3 5.5 
4 7.75 
5 10.6 
q ? 13.4 
| 9 18.5 
12 23.6 
13 26.2 
n 14 27.8 
| 
1 
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The SELECT BEST FIT key generates the display below: 


SELECT BEST FIT 


EQUATION 
n B RES ERROR 
Y = AX 
1.33339 9.24943 
Y= A + BK 
-8.16271 2.00973 0.24293 
Y = REEXPCEXID 
4.62083 8.12632 4.34318 
Y= 120 + BEO 
8. 17276 -0.91092 38. 26587 
Y - Re BX 
30.06560 -94.45703 11.38286 
Y = A+ BXLOGOD 
-12.45100 14.62717 2.68537 
Y = AKXHB RSN 
1.85781 1.83141 9.289809 
Y = K/CA + BEX) 
8. 54364 -0.09494 1.87836 


EQUATION Y = à + BKX HAS MAXIMUM R-SQUARE 


-SUUARE 


8.94563 


EQUATION Y = AXX HAS MINIMUM MAXIMUM ABSOLUTE RESIDUAL 


MAX DEUIATION 


9,632842 


8.71408 


eJ 
e 
cen 
uo 
Ga 
Gi 


The R-square column indicates that the Y A+BX has the maximum R-square. Y = AX has 


the smallest maximum deviation. 


For data in the restricted range, four equations fit better than Y = A+B log X. 


The equations Y = AX and Y = A+BX are plotted, and residuals are calculated. 
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Y = ARX 


À = 
1.993304209 


me 22. 3 23 53 


R-SQUARE = 
8.997109927539 


RES ERROR 
0.249482946299 


MAX CABSCRESIDUAL)) 
8.630478955007 


— AS 2». P, 3 
) 


en je | 


KES NGS FOR ERUS ToN : Y = AKX 


' VESTIMATE RESIDUAL 
n 3 5.5 5.981712627 -9.481712626996 
4 7,75 >.97561683599 -8.225616935994 
5 18.6 3.96352104493 6.620478955807 
7 13.4 13, 957329463 -0, 55732946299 
t 3 18.5 17,945137881 8.5549621138012 
| 12 23.6 23. 426850508 -8, 326850507982 
13 26.2 25.920754717 V.279245783419 
14 27.8 27.914653926 -0.114658325379 
1 
f Plot 50: Stat Vol 1 3 -77 


TEST AND REGRESSIONS 


SIMPLE REGRESSIONS || 


Y = A + B£X 


as 
-0.162?87722385 


B = 
2.00972629521 


R-SQUARE = 
8.997185818379 

RES ERROR 
8.242931736722 


MAXCABSCRESIDUAL)) 
8.714076246334 


ke ia a FOR EAM : Y= A + BRX 


YESTIMATE RESIDUAL l 
3 5.5 5.86647116325 -9.366471163246 L| 
4 7475 7.87619745846 -0, 126197458454 
5 19.6 9.80592375367 8.714076246334 
? 13.4 13.9053763441 -0.505376344081 
9 18.5 17.9248289345 0.575171965494 Vi 
12 23.6 23.954007820 1 -9,354007820137 
13 26.2 25.9637341153 0.236265884651 
14 27.8 27.9734604106 -0.173460418557? 
U 
Trimming the data to range 15 to 45 would probably be helpful in determining which (JÅ 


equation fits best in this range. 


U 
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ONE SAMPLE ANALYSIS 


INTRODUCTION 


DESCRIPTION 


The one sample analysis program will accept complete or censored samples, sort, edit, or 
trim the data, and print descriptive statistics. Six theoretical cumulative distribution func- 
tions considered by the program are Normal, Log Normal, Logistic, Weibull, Exponential, 
and Uniform. The program will calculate and print the least squares estimates of the para- 
meters for the distributions, and will plot their probability axes. Additionally, cumulative 
and relative frequency histograms may be drawn, as well as the empirical cumulative dis- 
tribution function for the data against the six theoretical cumulative distribution functions. 


Theoretical population percentiles, Kolmogorov-Smirnov and Cramer-VonMises goodness of 
fit tests, confidence intervals for the mean, and t-test of the sample mean may be calculated. 
Sample cumulative distribution functions may also be compared against theoretical cumu- 
lative distribution functions at each data point. 


HARDWARE REQUIREMENTS 


A Tektronix 4050-Series Graphic System with a minimum 16K memory is required for 
execution of this program. 


If the Data Communications Interface (Option 1) is installed, minimum acceptable memory 
configuration is 24K bytes. 


TAPE STRUCTURE 


Tape file 45 on the program tape contains overlay instructions required for program execution. 
When the program tape has been replaced by another tape, necessary reinsertion of the pro- 
gram tape will be indicated by the alarm bell and a displayed message. 


Magnetic tape file assignments for this program are specified in the table below: 


FILE NUMBER CONTENT OF THE FILE 


User definable key and program overlay operating system 


45 Data entry, data display, data plot, simple statistics, t-test and 
confidence intervals, data sort, data onto tape, data from tape, add 
data, delete data, change data, censor data, and trim data routines 
47 Data sort, histogram, cumulative histogram, and empirical CDF 
routines 
48 Data sort, quantiles, percentiles, and distribution parameters 
49 Data sort and probability plots 
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PROGRAM LIMITS 


For sample sizes greater than 100, programming modification is required. The variable H2 
establishes the dimension of the data file, and must be altered at two locations to accomodate 
larger sample sizes. The variable is defined on lines 1690 and 3920 of file 46, and increasing 
the dimension in these locations will allow larger samples to be entered. 


Variable H3 is assigned a value of 33, and thus allows data storage on the internal tape. 
Modification of the value of this variable makes data storage possible on other devices. 
Variable H3 is located on lines 3800 and 3880 of the store and load data routines contained 


in file 46. 


METHODS 


Density functions for the distributions which may be fit using the probability plotting 


methods in this program are: 


DISTRIBUTION 


Uniform 


Exponential 


Normal 


Lognormal 


Weibull 


Logistic 


4-2 


1 (log X — u)? 
——Ó—À oxo(- a 
V 210? 20? X>0 


B exp [ — (A+ BX)] 
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INTRODUCTION 


Data entered from the keyboard is stored in a one dimensional array in internal memory, and 
is stored to tape only if the STORE DATA key is pressed. Editing takes place in memory 
only, but edited data may be stored on tape with the STORE DATA key. 


The LIST DATA and PLOT DATA keys initiate action based on the data in memory, whether 
or not the data has been sorted. 


Probability plotting is accomplished with a probability axis which is constructed for each 
distribution selected. From the plot, parameters of the distribution are estimated, and a 
"best" fitting line is drawn. If the assumed distribution is correct, the plotted points will 
tend to fall in a straight line. 


If the plotted line appears not to fit censored or trimmed data, then the assumed model is 
not correct even though the data appears roughly linear. From the best fitting line, esti- 
mates of population parameters are found. Using these estimates, percentiles of the various 
distributions may be obtained. The sample cumulative distribution function and the assumed 
population distribution function are then compared. 


For a sample entered into memory, let ry and ra represent points censored on the left and 
right respectively. Points trimmed on the left and right are represented by r4 and r,. Let 
F(.) be the cumulative distribution function of themodel. LetV, = F^! (P;) where 


r +i-% 


PE 
Ntr, try 
i= rat, co NA 


The following are plotted for the distributions: 


DISTRIBUTION PLOT 


Uniform (V;,X;) 

Normal (V, X) 

Lognormal (Vi, log (X;)) 

Logistic (X;, log (P./(1 —P;))) 
Weibull (log X;, log (-log(1 —P;))) 
Exponential (—log(1 —P;), X.) 


A “best” line is then fit to the data points. 
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CENSORED DATA 


Censored data may be defined as a known number of observations with both known and un- 
known values. Censoring is a sample property, not a population characteristic, and may be ex- 
emplified by the following: An experimenter attempts to estimate the mean weight of a pop- 
ulation. 100 persons are to be weighed as the sample, but the scale available can only weigh 
objects in the range 100 to 200 pounds. 81 of the sample persons weigh between the range 
limits; 7 weigh less than 100; and 12 weigh more than 200. The results constitute censored 
data with a left censor point = 100, and right censor point = 200. We may say that 7 ob- 
servations are censored on the left, and 12 observations are censored on the right. 


LE E 


i. E 


E, 


TRIMMING DATA 

The trim data operation causes data points which have values outside the desired range to 
be eliminated from calculations and plotting. When the TRIM DATA key is pressed, the 
number of points to be ignored on the left (too small) and on the right (too large) must be 
specified. 


i te 


[= 


If the data has been sorted, the trimmed values will be ignored during the calculation of 

statistics and the t-test. If the data has not been sorted, trimming the data will have no effect V 
on statistics and t-test. The effect of trimmed data on probability plots and distribution 
parameters is to increase the number of censored data points. 


E VE 


REFERENCES 
Chernoff, H., Lieberman, G. J., “Use of Normal Probability Paper”, Journal of the American 
Statistical Association, Vol. 49, 1954. 


T ua 


Chernoff, H., Lieberman, G. J., “The Use of Generalized Probability Paper for Continuous 
Distribution”, Annals of Mathematical Statistics, Vol. 27, 1956, pp. 806—818. 


| 


Hahn, G. J., Shapiro, S. S., Statistical Models In Engineering, New York: John Wiley and 
Sons Inc., 1967. 


Sarhan, A. E., Greenberg, B. E., Contributions To Order Statistics, New York: John Wiley 
and Sons Inc., 1962. 
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VARIABLES 


Scratch Array 

Scratch Array 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

The array which stores the data sample of size 
N 

Scratch Array 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Scratch 

Stores trimmed or censored values 

Flag values (x) are: 1 7 left censor; 2 - right 
censor; 3 = left trimmed; 4 = right trimmed; and 
5 = data is sorted 

Indicates which file is to be appended into 
memory, and which routine of the file is to be 
executed 
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VARIABLE USAGE 


H2 Dimension of data array in the data entry 
routine and the data entry from tape routine 

H3 Device address for data storage 

H (x,y) Stores parameters for distributions 


x values are: 1 = Normal; 2 = Log normal; 
3 = Logistic; 4 = Weibull; 5 = Exponential: 
and 6 = Uniform. 

y values are: 0 = parameters undetermined; 
1 = parameters stored in an array 


N Number of elements in the data sample 
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OPERATING INSTRUCTIONS 


FUNCTION KEYS 


KEY 


A OU N > 


ONE SAMPLE ANALYSIS 


SHIFT KEYS 
ADD DELETE CHANGE 
POINT POINT 


CENSOR 
DATA 


QUANTILES | PERCENT 


KEY FUNCTION 


ENTER DATA allows entry of data from keyboard. 
LIST DATA prints entire data sample, including trimmed data. 
PLOT DATA plots data values. 


SIMPLE STATISTICS|prints variable N, total, mean, median, variance, standard 
deviation, minimum, maximum, range, skewness, kurtosis, 
3d moment, 4th moment, and the coefficient of variation 
for trimmed and censored sample. * 


T-TESTS calculates t-test of hypothetical mean; calculates con- 
CI fidence intervals for the mean for trimmed and censored 
sample.* 


SORT DATA sorts complete sample to ascending sequence in the same 
array in which it is stored. 


HISTO plots a histogram of the trimmed sample with a normal 
overlay if the normal parameters are calculated. * 


CUMM HISTO plots a cumulative histogram of the trimmed and censored 
data with a cumulative normal overlay if the normal para- 
meters are calculated.* 


* If the data has not been sorted, the complete sample is used. 
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Key 


Key Name 


Function 


9 QUANTILES calculates the sample cumulative distribution function 
(CDF), the distribution CDF, and the difference between 
the two for each data value; prints the Kolmogorav- 
Smirnov and the Cramer-VonMises tests; prints the dis- 
tribution parameters. 

10 PERCENT calculates theoretical deviate value for given cumulative 
probabilities. 

11 STORE DATA stores complete data sample to tape. 

12 LOAD DATA allows entry of data from tape. 

13 ADD DATA adds data to data previously entered. 

14 DELETE POINT removes incorrect value from sample. 

15 CHANGE POINT allows correction of data value. 

16 CENSOR DATA provides for the entry of the number of unknown values 
in the sample. 

17 TRIM DATA excludes out of range data values from calculations. 

18 PROB PLOT plots the trimmed and censored sample for selected 
distributions. 

19 DIST PARAM displays the distribution parameters of the trimmed and 
censored data. 

20 EMPIR CDF plots the empirical cumulative distribution function versus 


PROGRAM LOADING 


This program may be loaded automatically from the directory, or manually from the key- 


theoretical distribution. 


€ 


board. 


AUTOMATIC LOADING FROM THE DIRECTORY 


1. Insert the program tape. 


2. Press AUTO LOAD. The program directory will be displayed with the message ENTER 
THE PROGRAM NUMBER YOU WANT. 


3. Enter 27; press RETURN. 


4. When the I/O light goes out, execute the program. 
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OPERATING INSTRUCTIONS 


MANUAL LOADING FROM THE KEYBOARD 


1. Insert the program tape. 


2. Enter FIND 45; press RETURN. 
3. Enter OLD; press RETURN. 


4. When the 1/0 light goes out, execute the program. 


PROGRAM EXECUTION 


One Sample Analysis data entry and program execution are described in the comprehensive 
example which follows. 
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EXAMPLES I 
Measurements of a contaminant, expressed as a percentage, have been taken from a chemical L 
compound. The measurements have been tabulated below. 

OBSERVATION NUMBER CONTAMINANT PER CENT L 
1 .0503 
2 .0231 
3 .0689 i 
4 .0250 
5 .0120 
6 .0968 L 
7 .0305 
8 .0128 
9 .0363 || 
10 .0139 
11 .0402 
12 .0171 L 
During 4 other attempts, contaminant was present but in an amount too small to be mea- | 
cu 


sured. 


What is the likelihood of various amounts of contaminant appearing in future batches of the 
compound? 


Press ENTER DATA (key No. 1), to establish the data sample. 


One Sample Analysis Data Entry 


Enter data values, one by one 


23 
as 
ES. ES UI EL AE ES EG 


When the 12th observation has been entered, press LIST DATA (key No. 2) to verify accuracy 
of the entries. 


Es 
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EXAMPLES 


One Sample Analysis Data 


In this example, the list reveals the following errors: entry number 3 is entirely incorrect and 
should be deleted; entry number 8 was entered as .035 and must be changed to .0305; and 
observation number 9 (.0363) was omitted and must be added. 


DELETE POINT (key No. 14) will delete the entered value when RETURN is pressed. 


Delete Data 


Enter data values to be deleted, one by one 


To change .035 to .0305, press CHANGE DATA (key No. 15), enter the old and new values 
values as shown, and press RETURN. 


Change Data 


Enter the value to be deleted 
and the new value to be inserted 


old value, new value .0350,.0305 


old value, new value 


ADD DATA (key No. 13) will enable addition of the omitted data. 


Add Data 


Enter additional! data values, one by one 
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Again, the LIST DATA (key No. 2) is pressed to list the entries which are now correct. 


One Sample Analysis Data 


Press CENSOR DATA (key No. 16), establish the 4 unknown observations at the low end 
of the data sample, and press RETURN. 


Censor Data 


Left censored2=8 
Right censored=0 


enter new values 4,8 
censored routine complete 


SORT DATA (key No. 6) will cause the sample to be sorted to ascending sequence. A dis- 
played message acknowledges completion of the sort. 


Sort Data 
Data sort complete 
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PLOT DATA (key No. 3) will graphically display the data values. 


One Sample Analysis Data Plot 


8.104 


9.088 


8.072 


8.856 


0.04 


0.024 


8 2 4 6 8 18 12 


Probability plots will determine the distribution which best describes the data sample. 
PROBABILITY PLOTS (key No. 18) offers a choice of plots as shown on the displayed 
menu. 


Probability Plots 


Enter 
1 - Normal Plot 
2 - Log Normal 
3 - Logistic 
4 - Weibull 
% - Exponential 
6 - Uniform 


To illustrate this example, all 6 options are selected. 
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8.184 


8.872 


8.856 


8.024 


8.088 


Normal Plot 


8.81 1 20 46 60 86 


Censored: Lt= 4 Rt= 6 


Trimmed: Lt= 


99 99.93 


8 Rt: 0 
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Log Normal Plot 


6.91 
2 28 40 60 88 98 


Censored: Lt= 4 Rt= 6 Trimmed: Lt= 8 Rt= 8 
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EXAMPLES S 
I 
iy Logistic Plot | 
| 
88 
68 
48 LI 
20 A i 
WH 
! L 
8.008 8.024 0.64 8.056 8.072 0.088 8.104 
Censored: Ltz 4 Rt= 86 Trimmed: Lt= 0 Rt= 6 i 
i 
L 
i 
v 
L 
4 E 1 6 Plot 50: Stat Vol 1 
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ONE SAMPLE ANALYSIS 


EXAMPLES 
Weibull Plot 
99.9 
70 
38 
1 
8.81 
1.@E-4 
-4 -3 -2 
LNCA 
Censored: Lt= 4 Rt= 6 Trimmed: Lt= 8 Rt= 8 
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EXAMPLES V i 
U 
-— Exponential Plot | 
8. 888 L 
L 

8.872 E 
i 
8.856 
L 
0.04 x | 
Wi 
0.024 | 
8.008 i 


40 98 9? 


Censored: Lt= 4 Rt= 6 Trimmed: Lt= 8 Rt: 6 


[Ilo ES 


U 
| 
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t Uniform Plot 
! 8. 104 
0.088 
n 2.072 
n 8.056 
|! 
n 8.84 
MAN 
E o 8.824 
I 9.008 
28 48 68 88 198 
1 Censored: Ltz 4 Rt: 8 Trimmed: Lt= 9 Rt= @ 
The data is best fitted by the Weibull and Exponential plots, and the remainder of the an- 
7 alysis may be limited to those options. 
DISTRIBUTION PARAMETERS (key No. 19) will display the parameters used in plotting. 

f Distribution Parameters 
n Distribution Å B 

Normal 0.02147572490986 0.0740626208205 

Log Normal -1.706467259023 0.3977441020907 
Å Logistic -1.01754245687? 48.8932884929 

Weibull 41.7346112 1.085146163663 

/ N Exponential . 6.831854671E-4 0.9277119266234 
" EE Uni form -0.0265454545455 0.104583092909 
Censored! Lt= 4 Rt= 6 Trimmed: Lt= O Rt» 9 4-19 
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0.0968 


Kolmogorov-Smirnov Test 
Cramer-Von Mises Test 
Å 


Censored: 


Quanti les 


Enter 
1 - Normal Distribution 
2 - Log Normal 
3 - Logistic 
4 - Weibull 
5 - Exponential 
6 - Uni form 


Quantiles 


Sample CDF 


Lt= 


0.28125 
. 34375 


4 


Rt= 


Weibull 


0. 328920689569 
0. 347454594705 
0.372197719097 
9. 439452548855 


21211413181 
, 758764202284 
IO IET 
91 
«972 


Trimmed: Lt= 


8 


For the normal distribution A is the mean and B is the standard deviation. For the log 
normal A is the mean and B is the standard deviation of the logarithm of the data. 


The QUANTILES key (key No. 9) will cause a menu to be displayed, and the Weibull! and Ex- 
ponential options are selected. 


Difference 


-0.0476706885691 
-0.00270459470511 


8. 0340522809029 
0.0292974519452 


-8.01678602998644 


0.0156604642511 
B. 001446048615 


-0.00246141318088 


6. 0224857977156 
à. 009657 


Rt: 8 
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n EXAMPLES 
n Quantiles 
Data Sample CDF Exponential Difference 
1 8.012 8.28125 0.335269398116 -0.8540193981156 
- 0.0128 0.343? 8.3541847388913 -0.0164347889135 
0.0139 0.40625 8.379317728175 0.0269322718252 
0.8171 8.46875 8.447086877595 9.8217431224855 
ti 8.0231 8.53125 0.554662558996 -0.,0234125583962 
[d 825 8.59375 ð. 584172815532 ð. 80957718446763 
0.0385 0.65625 0.659028379467 -9.80277837946732 
0.8363 0.71875 A. 72341917763 -0.00466917762958 
n 8.0402 0.78125 ð. 759728519117 0.021521480883 
8.0583 0.84375 8.833114918432 0.01806358815879 
8.0689 8.90625 8,914705132115 -0.88845513211456 
a 0.8968 8.96875 8. 968834086242 -8.408624157E-5 
l Kolmogorov-Smirnov Test =0,0540193981156 
Cramer-Von Mises Test =z, 012486701682 
À z0. 000683185463094 
! B 20.0277119266234 


2 


Censored: Lt= 4 Rt: 8 Trimmed: Lt= 8 Rt: 8 


— 


When the PERCENTILES (key No. 10) is pressed, we again elect the Weibull and Exponential 
distributions from the selection offered. 


Percentiles 


Enter 


par ROLE een | 


NAS GIN = 
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Normal Distribution 


Log Normal 
Logistic 
Weibull 
Exponential 
Uni form 
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ND ND MD ND NO MD (1 Ma 
9 49 00 DOTS ( CO CO FO FV O 
49 C 02 1 & & C Q & NI O & nm 


Censored: Lt= 


AD % C UI & C» & C» OI MO O — 


M MO AO LO AD Nr OJ pm 
AØ AD €0 ^4 01 G5 O C AT FO e GO 


Censored: Lt= 


4 


Population Percentiles 


Weibull 


8. 0800403506901 604 
decur EE 


0. 0202962806277 
-9635774768018 


2 
FEN PEN 
f 


ó Trimmed: Lt= 


Population Percentiles 


Exponential 


0.080710911254924 
0.000961699622806 
0.8012430414058? 
0.00138479070022 
8.00210462147343 
0, 093667232034201 
0. 01989162927 
a. 064432254604: 
G. 08370069441 
B. 10290914221 


6 
6, 12830132: 
8 


Trimmed: Lts 


Ló] 
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n | EXAMPLES 
n EMPIRICAL CDF (key No. 20) provides option choices from the display shown here. 
n Empirical CDF 

Enter 
1 1 - CDF with Normal Overlay 

2 - Log Normal 

3 - Logistic 

4 - Weibull 
N 5 - Exponential 
| 6 - Uniform 

? - COF wzo Overlay 

Empirical CDF 
| i Weibul! Over lay 
MAN 
n 0.8 
n T 
l 9.4 
n 8.2 
| 
i l 
9.008 9.824 0,04 0.056 0.072 0.898 80.184 

f Censored: Lt= 4 Rt= B Trimmed: Lt 8 Rt= e 
n Plot 50: Stat Vol 1 4-23 
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Empirical CDF 
1 Exponential Overlay 


Lt= 4 Rt 


M 
us 
8 
0.008 8.024 8.84 6.656 8.672 8.888 8.184 | 
Censored: Lt= 4 Rt= 8 Trimmed: Lt= 8 Rt= 9 
Press SIMPLE STATISTICS (key No. 4) for the following display: i 
Simple Statistics l 
N 212. 
Total 20,4269 
Mean 28.035575 
Median 20,027?5 || 
Variance z0, 000661511136364 
Std Dev =0, 02571908587936 
Min 29.812 
Max =9.8968 |i 
Range =y, 3848 j 
Skewness =1,25647952264 
Kurtosis =3,69661731253 
3rd Moment 20,0800167619646563 | 
4th Moment =8,00606135705113638 
Coeff of Var 28, 722975651261 


Trimmed: Lt= 0 t= 8 
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ONE SAMPLE ANALYSIS 


EXAMPLES 


A menu is displayed when HISTO (key No. 7) is pressed, and information is entered as shown 
in the example. The histogram is displayed when RETURN is pressed at the completion of 
the entries. 


Histogram 


Minz0.012 
Maxs8. 3968 
Range=0.0848 


Enter 
1 to specify & of cells 
2 to specify cell width 
3 to specify offset and cell width 


Enter. offset and cell width .01,.81 


Histogram 
4 
2 
8 
90.01 0.03 6.65 0.07? 9.89 
Censored: Lt= 4 Rt= ô Trimmed: Lt= 8 Rt= 8 
Plot 50: Stat Voi 1 4 -25 


ONE SAMPLE ANALYSIS 


EXAMPLES 


CUMM HISTO (key No. 8) is used in the same manner. 


Cumulative Histogram 


Mins0.012 
Haxs8, 0968 
Range=0. 0848 
Enter 
1 to specify 4 of cells 
2 to specify cell width 
= to specify offset and cell width 


Enter offset and cel! width .61,..61 


Cumulative Histogram 


2.81 2.83 9.85 0.07 
Censored: Ltz 4 kt= 6 Trimmed: Lt= 
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ONE SAMPLE ANALYSIS 


EXAMPLES 


t-Test calculation is initiated by pressing T-TEST (key No. 5). In addition to the traditional 
confidence intervals of 90, 95, and 99, other confidence intervals may be specified. 


T-Test & Confidence Intervals 
Enter hypothesized mean .à4 


Do you vish any other confidence 
intervals besides 98, 95, & 99 (Y or ND Y 
How many add! conf. levels 1 

Enter conf. levels 

99.9 


Hypothesized mean =8,04 


TaTest 270. 595984984598 
Std error 20, 00742465369901 


Conf. Level Lower Limit Upper Limit 


98. 0.0222411345474 6. 8489888654526 
95.8 0. 0192334296767 8.0513162702222 
99.0 0.0125157619196 9,0526342 260804 
99.9 0. 00263567686797 Ø. 8655143231321 


To store the data sample to tape on a previously marked data file, press STORE DATA 
(key No. 11). The message displayed on the screen will prompt removal of the program 
tape, and request the data file number. 


To Tape 


Insert data tape and enter the tape file number 6 


Remove the program tape; insert the correct data tape; enter the file number of the data 
file. In this example, data file number 6 was used. 


When the data has been stored to the tape file, the completion message shown below will be 
displayed. 
Re-insert program tape. 


To tape complete 
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To continue program execution, re-insert the STATISTICS, VOLUME 1 program tape. 


To retrieve this data sample later, press LOAD DATA (key No. 12). The message shown be- 
low prompts removal of the program tape and requires the entry of the data file number. 


One Sample Analysis Data Entry From Tape 


Insert data tape and enter the tape file number 


When the data tape has been inserted, enter the file number and press RETURN. 


After all data has been entered, the following message will appear. 


RE-INSERT PROGRAM TAPE. 
FROM TAPE COMPLETE. 


When the program tape has been re-inserted, analysis may begin. 
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re SOFTWARE CHANGE INFORMATION 


^ TEKTRONDE! propucr PLOT 50 | CHANGE REFERENCE 432093 — 
l COMMER O MUR RR STATISTICS VOLUME 1 DATE __ 10-27-77 


CHANGE: DESCRIPTION 


EFF SN B060558 


062-1854-00 


REASON FOR CHANGE: 
1. 8K byte machine is unable to enter data. 16K byte machine gets a MEM FULL during 


calculations. 
2. The regression routine does not correctly trap negative or zero x and y values. 
DESCRIPTION OF CHANGE: 
1. Delete line 488 in File 36. Delete lines 640, 960, 1630, and 1720 from File 40. 
Change line 110 in File 40 to read W0=95+205*(MEMORY>8000). 


2. The following changes correct the problem: 


5 
7 
E 
a 


) 


a. Change line 1240 in File 42 to read 
IF B2+B3<0 THEN 1270 

b. Add line 505 to File 43 which reads 
IF B2+B3<0 THEN 710 


c. Add line 1145 to File 43 which reads 


IF B2+B3<0 THEN 1210 


"Es 


d. Change line 1110 in File 44 to read 
IF B2*B3«0 THEN 1140 

Add line 2305 to File 44 which reads 
IF B2*B3«0 THEN 2460 

f. Add line 2895 to File 44 which reads 


IF B2+B3<0 THEN 2960 
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